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Abstract. We consider the stability of periodic gravity free-surface water 
waves traveling downstream at a constant speed over a shear flow of finite 
depth. In case the free surface is fiat, a sharp criterion of linear instability 
is established for a general class of shear flows with inflection points and the 
maximal unstable wave number is found. Comparison to the rigid-wall set- 
ting testifies that the free surface has a destabilizing effect. For a class of 
unstable shear flows, the bifurcation of nontrivial periodic traveling waves is 
demonstrated at all wave numbers. We show the linear instability of small 
nontrivial waves that appear after bifurcation at an unstable wave number 
of the background shear flow. The proof uses a new formulation of the lin- 
earized water-wave problem and a perturbation argument. An example of the 
background shear flow of unstable small-amplitude periodic traveling waves 
is constructed for an arbitrary vorticity strength and for an arbitrary depth, 
illustrating that vorticity has a subtle influence on the stability of free-surface 
water waves. 



1. Introduction 

The water-wave problem in its simplest form concerns two-dimensional motion 
of an incompressible inviscid liquid with a free surface, acted on only by gravity. 
Suppose, for definiteness, that in the {x, y)-Cartesian coordinates gravity acts in 
the negative y-direction and that the liquid at time t occupies the region bounded 
from above by the free surface y = rj{t] x) and from below by the flat bottom y = 0. 
In the fluid region {(x,y) :0 < y < ri(t;x)}, the velocity held {u{t;x,y),v{t;x,y)) 
satisfies the incompressibility condition 

(1.1) d.^u + dyv^0 

and the Euler equation 



(1.2) 




- 9-, 



where P{t; x, y) is the pressure and g > denotes the gravitational constant of 
acceleration. The flow is allowed to have rotational motions and characterized by 
the vorticity to — Vx — Uy. The kinematic and dynamic boundary conditions at the 
free surface {y = r]{t; x)} 

(1.3) V ^ dtr] + udxT] and P = Fatm 

express, respectively, that the boundary moves with the velocity of the fluid par- 
ticles at the boundary and that the pressure at the surface equals the constant 
atmospheric pressure Patm- The impermeability condition at the flat bottom states 
that 



(1.4) 



1) = at {y = 0}. 
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It is a matter of common experience that waves which may be observed on the 
surface of the sea or on the river are approximately periodic and propagating of 
permanent form at a constant speed. In case lu = 0, namely in the irrotational 
setting, waves of this kind are referred to as Stokes waves, whose mathematical 
treatment was initiated by formal but far-reaching considerations of Stokes |53| 
himself. The existence theory of Stokes waves dates back to the construction due 
to Levi-Civita ^36] and Nekrasov 50J in the infinite-depth case and due to Struik 
[53] in the finite-depth case of small-amplitude waves, and it includes the global 
theory due to Krasovskii [35J and Keady and Norbury . Stokes waves of greatest 
height exist [58], [46] and are shown to have stagnation at wave crests fS*. A nice 
survey on the existence of Stokes waves include [59] and [HI Chapter 10, 11]. In the 
finite-depth case, it is recently shown by Constantin [16] that there arc no closed 
paths in the Stokes waves and each particle experiences a slight forward drift. 

While the irrotational assumption may serve as an approximation under certain 
circumstances and has been used in a majority of the existing research, surface 
water-waves typically carry vorticity, e.g. shear currents on a shallow channel and 
wind-drift boundary layers. Moreover, the governing equations for water waves 
allow for rotational steady motions. Gerstner [27] early in 1802 found an explicit 
formula for a family of periodic traveling waves on deep waters with a particular 
nonzero vorticity. An extensive existence theory of periodic traveling water waves 
with vorticity appeared in the construction due to Dubreil-Jacotin 25] of small- 
amplitude waves. Recently, for a general class of vorticity distributions, Constantin 
and Strauss [18] in the finite-depth case and Hur [32] in the infinite-depth case ac- 
complished the bifurcation analysis for periodic traveling waves of large amplitude. 
Partial results on the location of possible stagnation are found in [20, 60 . 

Waves of Stokes' kind is one of the few exact solutions of the free-surface water- 
wave problem, and as such it is important to understand the stability of these 
solutions. In this paper, we investigate the linear instability of periodic gravity 
water-waves with vorticity. 

The stability of water waves in case of zero vorticity has been under research 
much by means of numerical computations and formal analysis, especially in the 
works of Longuet-Higgins and his coworkers. Numerical studies of stability of Stokes 
waves under perturbations of the same period, namely the superharmonic pertur- 
bations, indicate that |48j . |55j instability sets in only when the wave amplitude is 
large enough to link with wave breaking ([3S]) and small- amplitude Stokes waves are 
found to be linearly stable under the same-period perturbations. The instability of 
large Stokes waves and solitary waves is recently proved by Lin ( 44] , [43]), under the 
assumption of no secondary bifurcation which is confirmed numerically. MacKay 
and Saffman [45] considered linear stability of small-amplitude Stokes waves by the 
general results of the Hamiltonian system. The Hamiltonian formulation in terms 
of the velocity potential, however, does not avail in the presence of vorticity ex- 
cept when the vorticity is constant [61] . The analysis of Benjamin and Feir [11] 
showed that there is a "sideband" instability for small Stokes waves, meaning that 
the perturbation has a different period than the steady wave. The Benjamin-Feir 
instability was made mathematically rigorous by Bridges and Mielke [12]. 

Analytical works on the stability of water waves with vorticity, on the other 
hand, are quite sparse. A recent contribution due to Constantin and Strauss [19] 
concerns two different kinds of formal stability of periodic traveling water-waves 
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with vorticity under perturbations of the same period. First, in case when the 
vorticity decreases with depth an energy-Casimir functional TL is constructed as 
a temporal invariant of the nonlinear water-wave problem, whose first variation 
gives the exact equations for steady waves il7i- In ilSj, the second variation of 7i 
is shown to be positive for some special perturbations and the water-wave system 
is called 7i-formally stable under these perturbations. The use of such energy- 
Casimir functional in studying stability of ideal fluids is pioneered by Arnold ^ 
for the fixed boundary case. The second approach of [12] uses another functional J , 
which is essentially the dual of 7i in the transformed variables but not an invariant. 
Its first variation gives the exact equations for steady waves in the transformed 
variables ([I7j), which served as the basis in [18] for the existence theory of traveling 
waves. The jZ-formal stability then means the positivity of its second variation. 
The "exchange of stability" theorem due to Crandall and Rabinowitz ^21, applies 
to conclude that the jZ-formal stability of the trivial solutions switches exactly at 
the bifurcation point and that steady waves along the curve of local bifurcation 
are J^-formally stable provided that both the depth and the vorticity strength are 
sufficiently small. 

The main results. As a preliminary step toward the stability and instability 
of nontrivial periodic waves, we examine the linear stability and instability of fiat- 
surface shear flows. The linear stability of shear flows in the rigid-wall setting is 
a classical problem, whose theories date back to the necessary condition for linear 
instabihty due to Rayleigh 52J. We refer to [38l EH EH ESI |39] and references 
therein for historic and recent results on this problem. In f39| , Lin obtained linear 
instability criteria for several classes of shear flows in a channel with rigid walls, and 
in this paper we generalize these to the free-surface setting. More specifically, our 
conclusions include: (1) The linear stability of shear flows with no inflection points 
(Theorem 16. 4|) . which generalizes Rayleigh's criterion in the rigid- wall setting [52] 
to the free-surface setting; (2) A sharp criterion of linear instability for a class of 
shear flows with one inflection value (Theorem l4.2p : and (3) A sufficient condition of 
linear instability for a class of shear flows with multiple inflection values (Theorem 
16. 1[) including any monotone flows. Our result testifies that free surface has a 
destabilizing effect compared to rigid walls. 

Our next step is to understand the local bifurcation of small-amplitude periodic 
traveling waves in the physical space. While our setting is similar to [Hj in that 
it hinges on the existence results of periodic waves in [TS] via the local bifurcation, 
the choice of the bifurcation parameter and the dependence of other parameters on 
the bifurcation parameter and free parameters in the description of the background 
shear fiow are different. In our setting, it is natural to consider that the shear profile 
and the channel depth are given and that the speed of wave propagation is chosen 
to ensure the local bifurcation. The relative flux and the vorticity-stream function 
relation are then computed. In contrast, in the bifurcation analysis |18j in the 
transformed variables, the wave speed, as well as the relative flux and the vorticity- 
stream function relation are held fixed. In turn, the shear profile and the channel 
depth vary along the bifurcation curve. Lemma 12.31 establishes the equivalence 
between the bifurcation equation (|2.7p (equivalently [TS] (3.8)]) in the transformed 
variables and the Rayleigh system (|2.9p - (|2.10p to obtain the bifurcation results for 
a large class of shear flows. In addition, our result helps to clarify the nature of 



4 



VERA MIKYOUNG HUR AND ZHIWU LIN 



the local bifurcation of periodic traveling water-waves that it does not necessarily 
involve the exchange of stability of trivial solutions (Remark I4.14|) . 

Our third step is to show under some technical assumptions that the linear 
instability of the background shear flow persists along the local curve of bifurcation 
of small-amplitude periodic traveling waves (Theorem l5.1|) . An example of such an 
unstable shear flow is 

U{y) = asmb{y - h/2) for ye[0,h], 

where h,b > satisfy hb ^ tt and a > is arbitrary (Remark l5.2p . In particular, by 
choosing a and h to be arbitrarily small, we can construct linearly unstable small 
periodic traveling water-waves with an arbitrarily small vorticity strength and an 
arbitrarily small channel depth. This indicates that the formal stability of the sec- 
ond kind in [19j (see discussions above) is quite different from the linear stability of 
the physical water wave problem. Our example also shows that adding an arbitrar- 
ily small vorticity to the water-wave system may affect the superharmonic stability 
of small-amplitude periodic irrotational waves in a water of arbitrary depth. Thus, 
it is important to take into account of the effects of vorticity in the study of the 
stability of water waves. 

Temporal invariants of the linearized water-wave problem are derived and their 
implications for the stability of the water-wave system are discussed (Section 3.3). 
In case when the vorticity-stream function relation is monotone, the energy func- 
tional d^Ti in [19j is indeed an invariant of the linearized water-wave problem. 
Other invariants are also derived. However, even with these additional invariants 
as constraints, the quadratic form d^Ti. is in general indefinite, indicating that that 
a steady (pure gravity) water-wave may be an energy saddle. A similar observation 
was made by Bona and Sachs [8] in the irrotational case. Therefore, a successful 
proof of the stability for the full water-wave problem would require to use the full 
equations instead of just a few invariants. 

Ideas of the proofs. Our approach in the proof of the linear instability of free- 
surface shear flows uses the Rayleigh system (I4.1|) - (I4.2|) . which is related to that 
in the rigid- wall setting [39] . The main difference from [39] lies in the complicated 
boundary condition (j4.2p on the free surface, which renders the analysis more in- 
volved. The instability property depends on the wave number, which is considered 
as a parameter. As in the rigid-wall setting [33], the key to a successful instability 
analysis is to locate the neutral limiting modes, which are neutrally stable solution 
of the Rayleigh system and contiguous to unstable modes. For certain classes of 
flows, neutral limiting modes in the free-surface setting are characterized by the 
inflection values. This together with the local bifurcation of unstable modes from 
each neutral limiting wave number gives a complete knowledge on the instability 
at all wave numbers. 

The instability analysis of small-amplitude nontrivial waves taken here is based 
on a new formulation which directly linearizes the Euler equation and the kinematic 
and dynamic boundary conditions on the free surface around a periodic traveling 
wave. Its growing-mode problem then is written as an operator equation for the 
stream function perturbation restricted on the steady free-surface. The mapping 
by the action-angle variables is employed to prove the continuity of the operator 
with respect to the amplitude parameter. In addition, in the action-angle variables, 
the equation of the particle trajectory takes a very simple form. The persistence of 
instability along the local curve of bifurcation is established by means of Steinberg's 
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eigenvalue perturbation theorem [51j . In addition, growing-mode solutions is proved 
to acquire regularity up to that of the steady profiles. 

This paper is organized as follows. Section 2 is the discussion on the local 
bifurcation of periodic traveling water-waves when a background shear flow in the 
physical space is given. Section 3 includes the formulation of the linearized periodic 
water-wave problem and the derivation of its invariants. Section 4 is devoted to the 
linear instability of shear flows with one inflection value, and subsequently. Section 
5 is to the linear instability of small-amplitude periodic waves over an unstable 
shear flow. Section 6 revisits the linear instability of shear flows for a more general 
class. 

2. Existence of small- amplitude periodic traveling water-waves 

We consider a traveling- wave solution of (|l.ip - (|1.4p . that is, a solution for which 
the velocity field, the wave profile and the pressure have space-time dependence 
{x — ct, y), where c > is the speed of wave propagation. With respect to a frame of 
reference moving with the speed c, the wave profile appears to be stationary and the 
flow is steady. The traveling-wave problem for (|l.ip - (|1.4p is further supplemented 
with the periodicity condition that the velocity field, the wave profile and the 
pressure are 27r/a-periodic in the a;- variable, where a > is the wave number. 

It is traditional in the traveling-wave problem to introduce the relative stream 
function ip[x^y) such that 

(2.1) %l)^ = -v, Tpy=U-C 

and ip{0, '7(0)) = 0. This reduces the traveling-wave problem for ()I.H) - ()I.4|) to a 
stationary elliptic boundary value problem 18, Section 2]: 

For a real parameter B and a function 7 G C^+'^([0, Ipol]), G (0, 1), find ri{x) 
and ip{x,y) which are 27r/Q;-periodic in the x-variable, ijjy{x,y) < in {{x,y) : < 
y < r]{x)J^ and 

(2.2a) -A'ip = j{ip) in 0<y<f]{x), 

(2.2b) il^^O on y^r]{x), 

(2.2c) |VV'|'+2gy = B on y^vi^), 

(2. 2d) il^^~Po on y = 0, 

where 

(2.3) Pa^ '4'vix,y)dy 

Jo 

is the relative total flu50. 

The vorticity function 7 gives the vorticity-stream function relation, that is, 
Lu — ^{ijj)- The assumption of no stagnation, i.e. tpy{x,y) < in the fluid region 
{{x,y) : < y < ri(x)}, guarantees that such a function is well-defined globally; 
See |18) . Furthermore, under this physically motivated stipulation, interchanging 
the roles of the y-coordinates and offers an alternative formulation to (|2.2|) in 



^In other words, there is no stagnation in the fluid region. Field observations |37) as well as 
laboratory experiments 1571 indicate that for wave patterns which are not near the spilling or 
breaking state, the speed of wave propagation is in general considerably larger than the horizontal 
velocity of any water particle. 

^PO < is independent of x. 
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a fixed strip, which serves as the basis of the existence theories in [25], [18], [32] , 
The nonhnear boundary condition (12.2b ) at the free surface y = r]{x) expresses 
Bernoulh's law. The steady hydrostatic pressure in the fluid region is given by 

(2.4) P{x,y)=B-^\Wi^{x,y)\^-gy- l{-p)dp. 



In this setting, a and B are considered as parameters whose values form part of 
the solution. The wave number a in the existence theory is independent of other 
physical parameters and hence is held fixed, while in the stability analysis in Section 
[5|it serves as parameter. The Bernoulli constant B measures the total mechanical 
energy of the flow and varies along a solution branch. 

2.1. The local bifurcation theorem in [18]. This subsection contains a sum- 
mary of the existence result in [TH] via the local bifurcation theorem of small- 
amplitude travelling- wave solutions to (j2.2p , provided that the total flux po and the 
vorticity-stream function relation 7 are given. 

A preliminary result for the local bifurcation is to find a curve of trivial solutions, 
which correspond to horizontal shear flows under a flat surface. As in [18|, Section 
3.1], let 

r(p) - f l{~p')dp'. r,„i„ = min r(p) sc: 0. 

Lemma 2.1 ([18 , Lemma 3.2). Given po < and e C^+^{[0, \po\]), (3 G (0, 1), 
for each e (— 2rniin,oo) the system i2.'2^) has a solution 



y{p) = 



dp' 



Po 



which corresponds to a parallel shear flow in the horizontal direction 



(2.5) u{t; X, y) = U{y- ^i) = c-^^x + 2r(p(y)) 

and v(t; x,y) = in the channel {{x, y) : < y < ft,(/i)}, where 

dp 



hifi) 



v//^ + 2r(p)' 



The hydrostatic pressure is P{y) — —gy for y € [0, /i(/i)]. Here, p{y) is the inverse 
of y = yip) o-nd determines the stream function ip(]j]^) — —p(y;^i); c > is 
arbitrary. 

In the statement of Theorem 12.21 below, instead of B the squared (relative) 
upstream flow speed /i = {U{h) — c)^ of a trivial shear flow (|2.5p serves as the 
bifurcation parameter. For each /i G (— 2rmin, 00) the Bernoulli constant B is 
determined uniquely in terms of /i by 

(2.6) B = n + 2g 



dp 

Po 



vV + 2r(p) 



The following theorem [TSJ Theorem 3.1] states the existence result of a one- 
parameter curve of small-amplitude periodic water-waves for a general class of vor- 
ticities and their properties, in a form convenient for our purposes. 
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Theorem 2.2 (Existence of small-amplitude periodic water-waves). Let the speed 
of wave propagation c > 0, thefluxpo < 0, the vorticity function -f £ C^^^([0, |po|]); 
(3 G (0, 1), and the wave number a > be given such that the system 

' {a^{iJ,)Mp)p = a^a{fi)M for p e (po, 0) 

(2.7) <^ M'/'Mp(0) = gM(0) 

M{po)=0 

admits a nontrivial solution for some /io G {-~2T -n-nn, oo) , where a(/i) = a{^;p) = 

v/M + 2r(p). 

Then, for e ^ sufficiently small there exists a one-parameter curve of steady 
solution-pair fi^ of i2.0j) and {r]^{x),ipeix,y)) of 112. 2\) such that r]^{x) and ipeix^y) 
are 2tt / a-periodic in the x-variable, ofC^^^ class, where (3 G (0, 1), and ijjty{x,y) < 
throughout the fluid region. 

At e = the solution corresponds to a trivial shear flow under a flat surface: 

(iO) The flat surface is given by rio{x) = /i(/io) =: ho and the velocity field is 

{ipOy{x,y),'-ipo^{x,y)) = {U{y) -c,0), 

where U (y) is determined in i2.5\) : 
(iiO) The pressure is given by the hydrostatic law Po{x, y) = —gy for y G [0, ha]. 

At each e > the corresponding nontrivial solution enjoys the following proper- 
ties: 

(ie) The bifurcation parameter has the asymptotic expansion 

/ie = /io + 0(e) as e ^ 

and the wave profile is given by 

ri^{x) — hf + a~^6je cos ax -\- 0{e'^) as e — > 0, 

where h^ — h{pL^) is given in Lemma \2.1\ and 5^ depends only on 7 and pq; 
The mean height satisfies 

= hci + 0{e) as e — > 0; 

Furthermore, the wave profile is of mean-zero; That is, 

(2.8) / {v,{x) - h,)dx = 0; 



(iie) The velocity field {il)^y{x,y),—il}^x{x,y)) in the steady fluid region {{x,y) : 
0<a;<27r/a,0<j/< //^(x)} is given by 

ipexix.y) = etjj^xiy) sinax -\- 0{e^), 

■>pey{x,y) ^ U{y) - c + etptyiy) cos ax + O(e^) 

as e —f 0, where ip^x o-nd ip^y are determined from the linear theory; 
(iiie) The hydrostatic pressure has the asymptotic expansion 

Pe{x,y) ^ -gy + 0{e) as e 0. 

The condition that the system (|2.7p admits a nontrivial solution for some /io G 
(— 2rmin,oo) is necessary and sufhcient for the local bifurcation [HI Section 3]. A 
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sufficient condition (^ISj) for the solvability of (12. 7|) and therefore the local bifurca- 
tion is 

/ (a^p ~ Po)'(2r(p) - 2r,nin)'/' + (2r(p) - 2r^i„)3/2) dp < gpl 
which is satisfied when po is sufficiently small. 

2.2. The bifurcation condition for a given shear flow. Our instability anal- 
yses in Section 4 and Section 5 are carried out in the physical space, where a 
shear-flow profile and the water depth are held fixed. The bifurcation analysis in 
the proof of Theorem 12.21 on the other hand, is carried out in the space of trans- 
formed variables, where the travel spped c, as well as the relative flux po and the 
vorticity-stream function relation 7 are held fixed, and the shear flow U (y) — c and 
the water depth h vary along the curve of local bifurcation. In this subsection, we 
study the local bifurcation in the physical space, with a given shear flow U (y) and 
the water depth h, which is relevant to the later instability analyses. The natural 
choice for parameters is the speed of wave propagation c > max U and the wave 
number k. 

Our first task is to relate the bifurcation equation (j2.7p in transformed variables 
with the Rayleigh system in the physical variables. 

Lemma 2.3. For the shear flow U{y) with y G [0, h] which is defined via Lemma 
\2.1\ the bifurcation equation 7| j is equivalent to the following Rayleigh equation 

(2.9) {U - c)(0" - k^(t>) - U"<j> = for ye (0, h) 
with the boundary conditions 

(2.10) + -r.d 

where {c,k) = (c, a), c > max [/, and (t>(jj) — {c~ U{y))M{p{y)). Here and in the 
sequel, the prime denotes the differentiation in the y-variable. 

Proof. Notice that c > max U. Indeed, 

a(M;p) - Vl^ + 2r(p) = ~iU{y{p)) - c) > 0. 

Since |H = -^'(y) = -(C/(y) - c), it follows that dp = 9,f^ = -jj^^dy. Let 
M{p{y)) = $(y), then (|2?7)) is written as 

((C/-c)2$')'-a2([/-c)2$ = for ye(0,/i), 

(U - cf<P'{h) = g^{h) and $(0) = 0. 

Let <j){y) = (c- U{y))^{y), and the above system becomes (|2.9p - (|2.10p . □ 

Remark 2.4. We illustrate how to construct downstream-traveling periodic waves of 
small-amplitude bifurcating from a fixed background shear-flow U{y) for y e [0, h]. 
First, one finds the parameter values {c,k) = (c, a) with c > maxt/ and a > 
such that the Rayleigh system (12. 9p - (|2.10p admits a nontrivial solution. The wave 
speed then determines the bifurcation parameter via — {U{h) — c)^, and the flux 
and the vorticity function are determined by 

(2.11) Po= {U{y)-c)dy and l{p) = U' {y{p)), 

Jo 
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respectively. By Lemma 12.31 the bifmcation equation (|2.7|) with /io, Po and 7 as 
above has a nontrivial sohition. Moreover, each shear flow U{h) — c for c > maxJJ 
corresponds to a trivial solution in Lemma |2. II Indeed, fi and c has a one-to-one 
correspondence via /i = {U{h) — c)^; The (relative) stream function defined as 

V'(y) = - / {U{h)~c)dy 



is monotone and its inverse y = y{—^) is well-defined. Then, Theorem 12.21 applies 
in the setting above to give a local curve of bifurcation of periodic waves. 

The lemma below obtains for a large class of shear flows the local bifurcation by 
showing that the Rayleigh system (I2.9p - (l2.10p has a nontrivial solution. 

Lemma 2.5. // 

(2.12) UeC^{[0,h]), U"{h)<0 and U{h) > U{y) for y ^ h, 

then for any wave number fc > there exists c{k) > U{h) — max[/ such that the 
system I12.9\) - I[2.10\) has a nontrivial solution (j) with (j) > in (0, h]. 

Proof. For c G {U{h), 00) and fc > 0, let <j)c be the solution of (|2.9p . or equivalently, 

(2.13) {{U -c)(j)'^-U'(j)c)' -k^{U -c)(j)e = for y G (0, /i) 

with (f>c{0) = and 0^(0) = 1. An integration of the above equation on the interval 
[0, h] yields that 

{U{h) - c)4>'Ah) - ([/(O) - c) - c^,{h)U' (h) - fcM ([/ - c)q^,dy = 



Note that the bifurcation condition (|2.10p is fulfilled if and only if the function 

(2.14) /(c) - c - U{0) + e {c~ U)4>cdy jTTr-.Mh) 

Jo c-U[h) 

has a zero at some c(fc) > U{h). It is easy to see that / is a continuous function of 
c for c > U{h). 

First, we claim that </>c(y) > for y € (0, /i]. Suppose, on the contrary, that 
(pcivo) — for some yo G (0, h]. Note that (|2.9p can be written as a Sturm-Liouville 
equation 

(2.15) - k"^, - 7^0c = 0. 

U ~ c 

Since 

[c-Uy + ^^{c-U) = Q for yG(0,/i), 

by Sturm's first comparison theorem the function c— U{y) must have a zero on the 
interval (0,yo)- A contradiction then proves the claim. 

Our goal is to show that /(c) > for c large enough and /(c) < as c — > [/ {h) + . 
Then by continuity, / vanishes at some c > U (h). First, when c — > 00 the sequence 
of solutions (j)c of (|2.9p . or equivalently (|2.15p converges in C^, that is, (f>c 4>oo- 
The limit 0oo satisfies the boundary value problem 

C - fc'<^oo =0 for y G (0, h) 

with 0oo(O) = and 0'oo(O) = 1. Therefore, (fioa is bounded, continuous, and positive 
on (0, h]. By the definition p.l4p . then it follows that /(c) ^ 00 as c ^ 00. 
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Next is to examine /(c) as c ^ U{h) + . Denote e = c — U{h) > 0. We claim 
that: 

(2.16) (t)c{h) ^ Ci > 0, for e > sufficiently small, 

where Ci > is independent of e. To see this, it is convenient to write (|2.13p as 
((c-C/)0^-(c-C/)'0,)' = fc2(c-C/)0c>O for ye{0,h), 

whence 

(2.17) (c-t/(y))0^(y)-(c-C/(y))V,(y) >c-f/(0) >0 for ye{0,h). 
Since 

(2.18) (c - C/)0^ - (c - C/)Vc = (c - [/)' 
it follows from ||TT7)) that 



U 



Mh)>ic-U{h)) I ,^ J' dy. 



,c-u J [c-uy 

and an integration of the above on [0, h] yields that 

c~U{0) 
lo {o-U{y)f 

Our assumption on U[y) asserts that Q ^U{h) — U{y) ^ P{h~y) for y e [h — 6, h] 
and S > sufficiently small, where /3 > is a constant. Thus, 

Mh) > e{c - Um f -— — i —dy 

Jh-5 {£ + (i[h~y)Y 

(l + Px) 

where Ci > is independent of e > 0. This proves the claim (|2.16p . 

To prove that /(c) < as c U{h)+, we consider the following two cases. 
Casel: max 0c (y) = 4'c{h). It follows that 



= (c - U{0)) I , ^ dx ^ Ci > 0, 




f(c)^c-UiO) + Ci\k' I ic-U)dy-—jj^j<0 

provided that c — U{h) = e > is small enough. 

Ca se 2: maxc/)^ = (ficiyc), where yc G (0, h). Since (A"(yc) and 0c(yc) > 0, by 
(|2T5l) we have U"{yc) > 0. On the other hand, U"{h) ^ and hence yc € [0,h- S] 
for some S > 0. Note that c — U{y) is bounded away from zero on the interval 
y G [0,h — S]. Since the coefficients of (|2.15p are uniformly bounded for c on 
y (z [0,h — d], the solution (pc is uniformly bounded on y e [0, /i — (5]. In particular, 
< 4>c{yc) ^ C2 independently for c. Therefore, 

/(c) ^ c - C/(0) + k^hC2 max(c - U[y)) - < 0, 

[OM] c — U yh) 

when £ = c — U (h) is small enough. This completes the proof. □ 



Lemma [2.31 and Remark 12.41 ensure the local bifurcation from a shear flow satis- 
fying p.l2p at any wave number fc > 0, as stated below. 
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Theorem 2.6. If U e C^{[0,h]), U"{h) < and U{h) > U{y) for y ^ h then 
for an arbitrary wavelength 27r/fc, where k > 0, there exist small- amplitude periodic 
waves bifurcating in the sense as in Theorem from the flat-surface shear flow 
U{y), where c{k) > maxC/. 

In the irrotational setting, i.e. [/ = 0, the parameter values c and k for which 
the Rayleigh system (|2.9|) - (|2.10|) is solvable give the dispersion relation (i.e. [22 j ) 

2 (7tanh(fc/i) 
" k ■ 

In the case of a nonzero background shear flow, such an explicit algebraic relation 
is in general unavailable. Still, the solvability of the Rayleigh problem (|2.9p - (|2.10p 
may be considered to give a generalized dispersion relation. Moreover, we have the 
following quantitative information about c(fc). 

Lemma 2.7. Given a shear flow U{y) in [0,ft.], let k and c{k) > maxt/ be such 
that i2.9\) ~ ^.10\) has a nontrivial solution. Then, 

(a) c{k) is bounded for k > 0; 

(b) If ki 7^ fc2 then c{ki) ^ cik-i); 

(c) In the long wave limit k —>■ 0+, the limit of the wave speed c(0) satisfies 
Burns condition 1151 

^^•'^^ ' {U-c{OW 9 

Proof. As in the proof of Lemma [^T5l the solvability of (|2.9p - (|2.10p is equivalent to 
the vanishing of the function / (c) defined by (|2.14p . 

(a) The proof of Lemma 12.51 implies that for each ^ > there exists Ca such 
that /(c) is positive when c > Ca, k < A. In interpretation, c(/c) ^ Ca for k < A. 
Thus, it suffices to show that /(c) > when c and k are large enough. Indeed, let 



dy 1 



c > max U be large enough so that 
solutions of 



u" 



^ 1 and let us denote by and 02 the 



+ (1 - k^) 01 = and (f)'^ + (-1 - fc2) 02 0, 

respectively, with 0^(0) — and 0^(0) — 1, where j — 1, 2. It is straightforward to 
see that 

(pi{y) = - ^^l ^ sinh(\/ fc^ - l)y and (l)2{y) = —^^=smh{\/k^~hT)y. 

Sturm's second comparison theorem [31] then asserts that the solution 0c, a- of (|2.15[) 
and 01 02 satisfy that 

01 ^ €,k ^ 02 
01 ^ 0c,fc 02' 

and thus 0i ^ 0c, fc ^02- It is then easy to see that /(c) > when k is big enough. 

(b) Suppose on the contrary that c(/ci) = c(/c2) = c for ki < k2. Let us denote by 
0c,fci and 0c,fe2 the corresponding nontrivial solutions of (|2.9p - (|2.10p . By Sturm's 
second comparison theorem IHT| , it follows that 

0c,fei(^) 0c,fe2(^) 

This contradicts the boundary condition (|2.10p . 
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(c) The Rayleigh equation (I2.13P for fc = implies that 



{c-U{y))ct>'M + U'{v)Uv)=m, 



where m is a constant. So by (|2.18p . 



and an integration of above from to ft, gives 



<l^c jh) 1 



c-U{h) "Vo (c-C/(y))2 



On the other hand, by (P?TU|) . 



m^{c- U{h))4>'^{h) + U'{h)(t>c{h) = 




A combination of above gives (|2.19p . 



□ 



The hmiting parameter value /i which corresponds to the limiting wave speed 
c(0) gives the lowest hydraulic head B defined in (|2.6p : see [18l Section 3] for detail. 
The limiting wave speed c(0) is the critical parameter near which small solitary 
waves of elevation exist [S5] , [33] . 

3. Linearization of the periodic gravity water-wave problem 

This section includes the detailed account of the linearization of the water- 
wave problem l|l.ip - (|1.4p around a periodic traveling wave which solves (|2.2p . The 
growing-mode problem is formulated as a set of operator equations. Invariants of 
the linearized problem are derived, and their implications in the stability of water 
waves are discussed. 

3.1. Derivation of the linearized problem of periodic water waves. A pe- 
riodic traveling-wave solution of (|2.2p is held fixed, and it serves as the undisturbed 
state about which the system (|l.ip - (|1.4p is linearized. The derivation is performed 
in the moving frame of references, in which the wave profile appears to be stationary 
and the flow is steady. Let us denote the undisturbed wave profile and (relative) 
stream function by rj^ix) and ipeixTu), respectively, which satisfy the system (j2.2p . 
The steady (relative) velocity field {ue{x,y) - c,Ve{x,y)) = {ipeyix,y), -^Jexix,!/)) 
is given by (|2.ip . and the hydrostatic pressure Pe{x, y) is determined by (|2.4p . Let 

Ve — {{x,y) : Q < X <2'K /a , Q < y < rje{x)} and Se — {{x^rje{x)) : Q < x < 2it / a} 

denote, respectively, the undisturbed fluid domain of one period and the steady 
wave profile. The steady vorticity ujf,{x^y) is given by lo,, ~ —Aipe — ^(ipe)- 

The linearization concerns a slightly-perturbed time-dependent solution of the 
nonlinear problem (ll.ip - p.4p near the steady state {rie{x),^e{x,y)). Let us denote 
the small perturbation of the wave profile, the velocity field and the pressure by 
r]e{x) + r]{t; x), (ue{x, y)-c + u{t; x, y), Ve{x, y) + v{t; x, y)) and Pe{x, y) -h P{t] x, y), 
respectively. We expand the nonlinear equations (|l.ip - (|1.4p around the steady state 
in the order of small perturbations and restrict the first-order terms to the steady 
domain and the boundary to obtain linearized equations for the the deviations 
ri(t; x), {u(t; x, y),v{t] x, y)), and P{t; x, y) in the wave profile, the velocity field and 
the pressure from those of the undisturbed state. 
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In the steady fluid domain Ve, the velocity deviation (m, v) satisfies the incom- 
pressibility condition 

(3.1) dxU + dyV = 

and the linearized Euler equation 



(3.2) 



dtu + {ue — c)dxU + UexU + Vf,dyU + v^yV ~ -^d^P 

dtV + (Ug - c)dxV + VexU + VedyV + V^yV = -dyP, 



where P is the pressure deviation. Equation (|3.ip allows us to introduce the stream 
function ip{t] x, y) for the velocity deviation {u{t] x, y), v{t; x, y)): 

dx''P = ~v and dytp = u. 

Let us denote by a;(t; x, y) the deviation of vorticity. By definition, lo = — Ai/j. The 
linearized vorticity equation is 

(3.3) dtUJ + {ipeydxUJ - IpexdyU) + {uJexOylp - LUeyOxlp) = 0. 

Since tUe = 7('0e), the last term can be written as {il^e)i''Peydx'4' ~ '^exdyt^). 

The linearized kinematic and dynamic boundary conditions restricted to the 
steady free boundary Se are 

(3.4) V + VeyTj dt-q + (Ue " CjOxTj + {u + U^yr])'qex 

and 

P + Peyr^ = 0, 

respectively. In terms of the stream function, (|3.4p is further written as 

= 9t?7 + ifjeydxTI + [dxi^ + Vexdylp) + (?/; 

exy + VexVeyyjV 

= dfll + Ipeydxfl + drijj + C^rV-'ea?? 
= dtl] + drilpeyV) + ^rlp, 

where 

drf = dxf + riexdyf 

denotes the tangential derivative of a function / defined on the curve {y = rie{x)}. 
Alternatively, drfix) — dxf{x,rie{x)). The bottom boundary condition of the 
linearized motion is 

dxil^ = on {y 0}. 

Our next task is to examine the time-evolution of tjj on the steady free surface 
Se- This links the tangential derivative of the pressure deviation P on the steady 
free surface Se with -0 and rj on Se- For a function / defined on Se, let us denote 

by 

dnf = dyf - rjexdxf 

the normal derivative of / on the curve {y = ?7e(a;)}. 

Lemma 3.1. On the steady free surface Se, the normal derivative dnip satisfies 

dtdnll^ + 5r(-0ey9„V) + f^^^V' + drP = 0, 

where fl = uje{x,rie{x)) is the (constant) value of steady vorticity on Se- 
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Proof. The linearized Euler equation ()3.2p can be rewritten in the form 

/ u \ I V 

(3.5) ~VP = dt{ ^ 1+V (K - C)U + VeV) -UJei^ J - I _ ^) 

by linearizing the nonlinear convection term according to the identity 

u \ /I / 2 , 2^^ f V 



Taking dot product of p.Sp with the vector (1, rjex) and restricting the result to 
Se, we have 

-drP = dt (u + Vriex) + dr {{Ue ~ C)u + VeV) - {v - Ur/ex) 

= dt (t/jy - rjexi'x) + dr {{Ue - C)%j}y - (Ug - c)r]exlj-'x) " ^ {-^x " Vexi'y) 
^ dtdnil + dr{{Ue-C)dn'4^)+^dr1p, 

where in the above derivation we use the steady kinematic equation 

= (Ue - d)Vex, on 5e. 

□ 

In summary, there results in the linearized water-wave problem: 

(3.6a) dfUJ + {TpeydxUJ - TpexdyUj) ^ l'ilpe){i>eydxll^ - i>exdy^) in T^e, 

(3.6b) dtT] + drii^eyV) + dri^ = On Se; 

(3.6c) P + Peyr]^0 on Se; 

(3.6d) dtdn'ip + drii^eydnll') + ^dri^ + drP ^ OU Se] 

(3.6e) dxiJ = on {y = 0}. 

Note that the above linearized system may be viewed as one for '>p{t;x,y) and 
r](t;x). Indeed, P{t;x,r]e(x)) is determined through p.6cp in terms of r]{t]x) and 
other physical quantities are similarly determined in terms of 'ip{t] x, y) and ri{t; x). 

3.2. The growing-mode problem. A growing mode refers to an exponentially 
growing solution to the linearized water-wave problem ()3.6I) of the form 

{in{t;x),^{t;x,y)) = ie^'r,{x),e^'ib{x,y)) 

and P{t;x,rie{x)) — e^* P (x , r]e{x)) with Re A > 0. For such a solution, the lin- 
earized vorticity equation (j3.6ap further reduces to 

(3.7) XUJ + (ipeydx^^ - Tpexdy^) - 7'{i^e)ii^eydxi^ - ^exdylp) ^ in , 

where lu = —Aip. Let {Xe{s;x,y),Ye{s;x,y)) be the particle trajectory of the 
steady flow 

iXe=^JeyiXe,Ye) 
\Ye = -^Jex{Xe,Ye) 



(3.8) 



with the initial position (Xe(0), Ye(0)) = {x,y)- Here, the dot above a variable 
denotes the differentiation in the s-variable. Integration of (|3.7p along the particle 
trajectory (Xe{s;x,y),Ye{s,x,y)) for s € (— oo,0) yields [41] Lemma 3.1] 

(3.9a) Ai^ + 7'(V'e)i^-7'(V'e) / Xe^'ij{Xe{s),Ye{s))ds = inVe. 
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For a growing mode the boundary conditions (I3.6b[) . (I3.6cp and (|3.6d[) on Se 

become 

(3.9b) \7]{x) + -^{ipeyix,r]eix))rj{x)) = --^7/>(a;,?7e(x)), 

(3.9c) P{x,7Je{x)) + Pey{x,7]e{x))r]{x) = 0, 

(3.9d) Xipnix) + -^{tpeyix,r]e{x))i)nix)) = - ^P{x , r]e{x)) ~ 17-^?/>(x, T^g (z) ) . 

The kinematic boundary condition p.6ep at the flat bottom {y = 0} implies 
that 0) is a constant. Observe that (|3.9aP - (j3.9dp remain unchanged by adding 
a constant to ip. So we can impose the bottom boundary condition by 

(3.9e) V(a;,0) = 0. 

In summary, the growing-mode problem for periodic traveling water-waves is to 
find a nontrivial solution of (|3.9ap - p.9ep with Re A > 0. 

3.3. Invariants of the linearized water-wave problem. In this subsection, the 
invariants of the linearized water-wave problem (|3.6ap - (|3.6eP are derived, and their 
implications in the stability of water waves are discussed. 
With the introduction of the Poisson bracket, defined as 

[/l, /2] = d^fidyf2 - dyfida:f2, 

the linearized vorticity equation p. Gap is further written as 

(3.10) 5tt^- [V'e,w]+7'(^e)[V'e>] =0 m2?e, 

where uj = —Aip. Recall that 

drf = dxf + rjexdyf and 9„/ = dyf - rjexdxf, 

where / is a function defined on Se- 

For the case when the vorticity-stream relation is monotone, we show that the 
following energy functional is an invariant of the linearized system. 

Lemma 3.2. Provided that either ^'{p) < or "/'{p) > on [0, \po\], then for any 
solution (rjjtp) of the linearized water-wave problem I13.6\) . we have ■^£{r],^p) — 0, 
where 



(3.11) 



f(?7,V')-//^ ^IVipfdydx- JJ^ ^j'iiPe)-'\co\^dydx 



- / ^Pey\ri\'^dx + Re / tpeydni^rj* dx - I ^ipeylvl'^dx. 

In the irrotational setting, i.e. 7 = 0, the invariant functional becomes 

(3.12) £(/?,?/')= // ilV^Apdydx- / \Pey\ri\'dx + Ke [ i^eydni^V* dx. 
J JVc Js^ 

Proof. By integration by parts, 

Jtjj \\^i^?dydx^Re jj dt{Vi^) ■ dydx 



= Re dtw>p*dydx+ / dt{dnip)ip* dx + / dtdy-tjjip* dx 

J Jv^ Js^ J{y=0} 

:= (/) + (//) + (///). 
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A substitution of dtOJ by the linearized vorticity equation (|3.10|) yields that 
(/)=Re JJ^ 4^] + [A, u:])rdydx 

= Re JJ^ {-^[iJen'ii'em^] + [i^e,u;r]-[i^e,r]to)dxdy. 
The second equality in the above uses that 

i[V'e,7'('^e)|V'P] = i[V'e,7'(V'e)]|^P + (V'e) ( [V-e , ^] V"* + [V'e,^*]^') 
= 7'(^e)Re[V'e,V']V'*, 

which follows since [tpe, g{ilJe)] = for any g. With the observation that 
J J [ipe, f]dydx = J J V • {{tpey, -tpex)f)dydx 

(V'ey, -^Jex) ■ {-Vex, i)fdx + [ f-ipexdx = 

for any function /, the integral (/) further reduces to 

(/) =Re // -[i>e,^]*u;dydx. 



v=o 



A simple substitution of [ipe,tp] by the linearized vorticity equation (jS.lOp then 
yields that 

(/) = Re // -f'iyje)-HdtCJ* - [i^e,uj]*)cjdydx 



(3.13) = Re {^j'{^/je)-'dt\iu\' - i[V'e,7'(V'e)-VP])rfydx 



d 

di 



^j'i^P,)-'\L,\^dydx. 



This uses that 



i[V'e,7'(V'e)-VP] = 7'(^e)"'Re[V'e,H*'^- 

Next is to examine the surface integral (//). Simple substitutions by the lin- 
earized boundary conditions (|3.6bp . (|3.6cp and (j3.6dp into (//) and an integration 
by parts yield that 

(//) = Re / driPeyV - ^eydn^P " dx 



Re / {PeyT] - 1peydn'ip){-dr'ip*)dx 

Re / (Pej,ry - 1peydn'^){dtVl* + drilpeyV*))dx 

/ ^Pey\ri\'^dx - Re / ipeydntpdt'ri*dx + Re / dr{P + tpeydnip)ipeyV*dx. 



d 

dt 



The second equality uses that 

"d-r = 

2 J c dx 



Re I {dr^)ilj* dx = \ I -^\ip{x,r]e{x))\^ dx = <d. 
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More generally, Re {dTf)f*dx = for any function / defined on Se- With 
another simple substitution by the boundary condition (|3.6dp . the last term in the 
computation of (//) is written as 



Re / dr{P + 1peydnij)4'eyV*dx 

= —Re / {dtdni^ + ^dr^)'lJjeyTI* dx 



= -Re / {dt{dnil^)^ey'n* dx + riRe / {dtrj + dr{^pey'n))ipey'n*)dx 

= -Re J ^eydtidni')r]*dx + Q-^ J ^^eyM'^dx. 
The last equality uses that Re dr{ipeyV){'4'eyV)*dx = 0. Therefore, 

(3.14) (//) = A£ ip,,y|,y|2da;_Rc^y" jP,ydn^V*dx + n^^ J \ijey\V?dx. 

Finally, it is straightforward to see that 

(///) ^'4,*[x,Q) I dt {dy^p)dx = 0. 
J{v=o} 

This, together with (|3.13p and (I3.14p proves that £{ri,^) is an invariant. In the 
irrotational setting, i.e. 7 = 0, the area integral (/) is zero, i7 = 0, and the other 
terms remain the same. This completes the proof. □ 

Remark 3.3. Our energy functional £ agrees with the second variation d^Ti. of the 
energy-Casimir functional in [19 . Recall that the hydrostatic pressure of the steady 
solution is given as in (12. 4p by 



ftpe{x,y) 

Pe{x,y) = B-^\VM^.y)\^ -9y + I i{-p)dp, 



where B is the Bernoulli constant. Differentiation of the above and restriction on 
Se then yield that 

-Pey - il^pey = ^Sj^lVV'eP +9, 

where fl = jiip = 0) = ujeix^rjeix)). Thus, when 7 is monotone, it follows that 



2f(77,V) = JJ^ l^^Pl^dydx- jj^ j\iJe)-^\u;\^dydx 

+ [ {hdy\y^e\^ +9)\V?dx + 2Re [ i^eydn^l^* dx. 

This is exactly the expression of the second variation d^Ti in [19] of the following 
enery-Casimir functional (in our notations) 

(3.15) n{v,tP)^ ff + gy ^ B - F{,u)] dydx, 



around the steady state (r]e,ipe)- Here, = {{x,y) : < y < r]{t]x)} and 
{F')~^ = 7. The quadratic form d^TC is used in [19] to study a formal stability. 

Our next invariant is the linearized horizontal momentum. The result is free 
from restrictions on 7. 
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^ / (V' + ^eyV)dx = 



Lemma 3.4. For any solution {-q, of the linearized problem of 113. 6\} . the identity 

d_ 
dt 

holds true. 

Proof. We integrate over the steady fluid region of tlie linearized equation for 
the horizontal velocity 

and apply the divergence theorem to arrive at 

(3.16) ^ J J dyTpdydx + j {dyi),-dxi))-{~riex,l)ipeydx^-JJ P^dydx. 
It is straightforward to see that 

dyipdydx — ip dx. 
In view of (j3.6bp . the second term on the left hand side of (|3.16p is written as 

{dr1p)lpeydx ^ / {dtrj + dr{lpeyTl))^eydx 

~ II J '^'^V'ndx- j tpeydr{iJey)'r]dx. 

With the use of Stokes' theorem and the dynamic boundary condition p.ip . the 
right side of p.l6p becomes 



Pxdydx = I P7]exdx = / Pey-qrjexdx. 

On the other hand, the steady Euler equation restricted to the steady free-surface 
Se yields that 

Pex ^ey'4^exy '4^ex'4^eyy — '4^ey{'4^exy "t" ^ea;^eyy) — ^ey^^r'^ey 

Since P^x + PeyVex = on 5e, a simple substitution then proves the assertion. □ 

Next, the integration of (|3.6bp on Se shows that rj dx is an invariant. Finally, 
multiplication on the linearized vorticity equation (|3.10p by ^ and then integration 
yield that 

j j uj£,dydx 

is an invariant, for any function 

£, e ker(V'ey9:r " ^^exdy) C L^iVe). 

We summarize our results. 
Proposition 3.5. The linearized problem US. 6\ ) has the energy invariant: 

f(r/,V)= JJ^ ^\Vi,\^dydx- JJ^ iy{A)-^\uj\^dydx 

+ I \{dy(\\V^,\'')+g)\rj\'dx + ^e f i;,yd^^r,*dx 
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when 7 is monotone and 

when 7 = (irrotational) . In addition, \3. 0jl has the following invariants: 



■>pey'n)dx, 



m{r], ijj) 



rjdx 



ui^dydx, for any f S ker('0oyC^x — 4'cxdy)- 



The nonlinear water-wave problem (|l.ip - p.4p has the following invariants [THl 
Section 2]: 



V{t) 



( i (u^ + ) + 5y) dydx (energy) , 



971 



m 



u dydx 



v{t) 



dydx 



V{t) 



f{uj)dydx 



(horizontal momentum), 



(mass). 



(Casimir invariant). 



V{t) 



where 

V{t) = {{x,y) -.0 <x < 27r/a , < y < r^it; x)} 

is the fluid domain at time t of one wave length, and the function / is arbitrary 
such that the integral ^ exists. 

The invariants of the linearized problem f , M., m, and J- in Proposition 13.51 
can be obtained by expanding the invariants of the nonlinear problem £, 971, m 
and 5^, respectively, around the steady state {rje{x),'ipe{x,y)). The quadratic form 
£{r],tp) is the second variation of the energy functional Ti defined in (|3.15p (see 
Remark 13.31) . which is a combination of €, 971, m and ^. The invariants A4 and m 
of the linear problem are the first variations of 971 and m, respectively. We note 
that is the first variation of 5", since by the assumptions of no stagnation and the 
monotonicity of 7 it follows that 

keT{'ijjeydx - ^exdy) = {/(V'e) : / IS arbitrary} = {f{i^e) ■ f is arbitrary}. 

We now make some comments on the implications of these invariants on the 
stability of water waves. A traditional approach of proving stability for conservative 
systems is the energy method, for which one tries to show that a steady state is 
an energy minimizer under the constraints of other invariants such as momentum, 
mass, etc. This method has been widely used in the stability analysis of various 
approximate equations such as the KdV equation [Q, [lOj and water waves with a 
large surface tension [17]. However, steady waves of the full pure gravity water- 
wave problem in general are not (constrainted) energy minimizers. This was first 
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observed (p|) in the irrotational case. By our discussions below, even with the 
favorable vorticity, the situation remains the same. Note that if a steady gravity 
water-wave is an energy minimizer under the constraints of fixed OJl, m and 3^, 
then at the linearized level the second variation £ should be positive under the 
constraints that the variations M.^ m and J- are zero. 

In the irrotational setting, the first two terms in the expression (|3.12p of £{rj^il)) 
yield a positive norm 

(3.17) // \\7tljfdydx+ I 

(Indeed, since APg = —2il)l^y — V'La; ~ ^ and Pey(a;,0) = —g by the 

maximum principle Pg attains its minimum at the free surface Se- Furthermore, 
since Pe takes a constant on Se by the Hopf lemma Pgj, < on iSg.) The last term 
Re tpeydn4'V*dx of the right side of (|3.12|) . however, does not have a definite 
sign and contains a 1/2-higher derivative than that of p.l7p . Thus, it cannot be 
bounded by the norm (|3.17|) . even with the constraints M = m = 0. Consequently, 
the quadratic form £{ri,il)) is highly indefinite unless ij^ey = 0, that is, for a trivial 
flow. In other words, steady gravity water waves in the irrotational case are in 
general not (constrainted) energy minimizers. Indeed, in [29], [13j . steady water- 
waves were constructed as energy saddles by variational methods. 

With a nonzero vorticity, a control of the mixed-type term in the energy func- 
tional by other positive terms fails for the same reason as in the irrotational setting. 
Let us consider the case of a favorable vorticity with 7' < 0. Under some additional 
assumptions, it is shown in [19] that the first three terms in the expression p.ll|) 
of £(77, ip) gives a positive norm 

(3.18) // (|VV'|2-H|c^ndydx+ / Ir^pdx. 

However, due to the lack of control of the boundary value of on Se, one could 
not use the elliptic regularity to get a better control for i/),such as |iV'llH2. Thus, 
the mixed-type term Re Jg tpeydnip'r]*dx is still not controllable by ()3.18|) . In [19], 
several class of perturbations are introduced to make this mixed term controllable 
by p.l8p . and therefore get the positivity off (77, tp) and a formal stability for these 
special perturbations. However, it is difficult to see that these special classes are 
invariant during the evolution of the water-wave problem. So, it remains unclear 
how to pass from such formal stability to the genuine stability, even for these 
special perturbations. It is not hard to see that by adding other invariants as 
constraints, one can relax the assumptions to get the positive term (|3.18|) such as in 
[5Tj for the fixed boundary case, but the mixed-type term remains uncontrollable. In 
conclusion, the quadratic form £{r], ip) is in general indefinite also in the rotational 
setting, and steady water waves with vorticity are expected to be energy saddles. 

The above discussions of steady water waves as energy saddles do not imply 
that steady water-waves are necessarily unstable. Indeed, as mentioned in the 
Introduction, the small Stokes waves are believed to be stable [IS], [5S] under 
perturbations of the same period. For the rotational case, under the assumption 
of a monotone 7, the corresponding trivial solutions with shear flows defined in 
Lemma [2T] can not have an inflection point, since j'{ipe{y)) = —U"{y)/U{y) 7^ 
and [/ < by the no-stagnation assumption. Thus by Theorem 16. 4( such shear 
flows are linearly stable to perturbations of any period. Since small-amplitude 
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waves with any monotone vorticity relation 7 bifurcate from these strongly stable 
shear flows, they are likely to be stable. But, a successful stability analysis of any 
nontrivial steady gravity water waves would require to use the full set of equations 
instead of a few invariants. 

4. Linear instability of shear flows with free surface 

This section is devoted to the study of the linear instability of a free-surface 
shear flow {U{y),0) with y E [0,h], which is a steady solution of the water-wave 
problem (|l.ip - (|1.4p with P{x, y) — —gy. In the traveling frame of reference with the 
speed c > max U, this may be recognized as a trivial solution of the traveling-wave 
problem ((^ : 

rjeix) = h and i^peyix,y), ~^JJexix.,y)) ^ {U{y) - c,0). 

Throughout this section, wc write U for U — cto simplify our notations. 
We seek for normal mode solutions of the form 

rjit; x) = r7,e^"(--=*) , i>{t-,x,y) = ,/.(y)e^"(--^*) 

and P{t\x,y) = P(y)e*"(^~^*) to the linearized water-wave problem p.6p . Here, 
a > is the wave number and c is the complex phase speed. It is equivalent to find 
solutions to the growing-mode problem (|3.9p of the form A = —iac and 

77(2;) -77^6^"", =0(y)e^"- 

and P{x, rie{x)) = P^e*""'. Note that RcA > if and only if Imc > 0. 

Since 7'(^e(j/)) = u:,y{y) /i^,y{y) = -U"{y)/U[y) and [X,{s),Y,{s)) = {x + 
U{y)s, y), the linearized vorticity equation p.9ap translates into the Rayleigh equa- 
tion 

(4.1) ([/- c) ((/)"- - = for ?;e(0,/i). 

Here and elsewhere the prime denote the differentiation in the y-variable. The 
boundary conditions (j3.9bp . (|3.9cp and (|3.9dD on the free surface are simplified to 
be 

ic-U{h))7jh^Hh), Ph-9Vh = 

and 

ic-U{h))q>'{h) = Ph + U'{h)cl){h), 
respectively. Eliminating ijh and Ph from the above yields that 

{U{h) - c)'cf>'{h) - (5 + U'{h){U{h) - c)) 

The bottom boundary condition p.9ep becomes (/)(0) = 0. So the linear instability 
is reduced to study the Rayleigh equation (j4.ip with the boundary condition 

/ (Uih) - c)^'{h) = (g + U'{h){U{h) - c)) m 
10(O) = O. 

A shear profile U is said to be linearly unstable if there exists a nontrivial solution 
of gT])-g21) with Imo 0. 

The Rayleigh system (|4.1l) - (|4.2p is alternatively derived in 62J by linearizing 
directly the water-wave problem (|l.ip - (|1.4p around (C/(y),0) in y € [0,/i]. Note 
that for a real number c > max[/ the Rayleigh system (|4.ip - (|4.2p becomes the 
bifurcation equation (p?9 l) - ([2T0| . 
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In case of rigid walls at y = h and y = 0, that is, the Dirichlet boundary 
conditions (l){h) = = (/){0) in place of (14. 2p . the instability of a shear flow is a 
classical problem, which has been under extensive research since Lord Rayleigh 
[52] . Recently, by a novel analysis of neutral modes, Lin [39] established a sharp 
criterion for linear instability in the rigid-wall setting for a general class of shear 
flows. Our objective in this section is to obtain analogous results in the free-surface 
setting. 

Below is the definition of the class of shear flows studied in this section. By an 
inflection value we mean the value of U at an inflection point. 

Definition 4.1. A function U £ C^([0, /i]) is said to be in the class /C+ if U has 
an unique inflection value Us and 

is bounded and positive on [0,h]. 

An example of class K.^ flows is U{y) = sin my, for which K{y) = m^. 
For U (y) in the class IC~^ , consider the Sturm-Liouville equation 

(4.4) (/)" - + K{y)(l3 = for y G (0, h) 

with the boundary conditions 

U'ih) = griUsmh) if U{h)^Us 
U(M=0 if U{h)^Us, 



(4.5) 



(4.6) (/)(0) = 0. 

Here, 



9 , U'ih) 



The following theorem gives a sharp instability criterion for shear flows in class 
/C+, for the free surface case. 

Theorem 4.2 (Linear instability of free-surface shear flows in /C"^). Consider a flow 
U (y) in class JC^ . Denote by —a^^^ the lowest eigenvalue of ~ K (y) with 
the boundary condition ^4-5^ ~ (T^ : which is assumed to be negative. Then for each 
a € (0, ttniax) , there exists an unstable solution-triple {(f), a, c) (with Imc > 0) of 
l{4.1\ )- (4^ - The interval of unstable wave numbers (0, amax) is maximal in the 
sense that the flow is linearly stable if either the operator ~ K{y) on y £ (0, h) 
with |^.5[)-(|7^ is nonnegative or a ^ ctmax- 



^From the usual variational consideration, the lowest eigenvalue — a^jj^^ 
acterized as 

, , 2 . , J,\\<l>'\'-K{y)mdy + gr{Usmh)\' 
(4.8a) -amax= mf — , , 

0(O)=O •'^ 

in case U{h) ^ Us, and 

4-8b) - "max = inf — rh, nn , 

0(O)=O=0(/.) -^0 ^ 
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in case U{h) — Us- 

4.1. Neutral limiting modes. The proof of Theorem 14.21 makes use of neutral 
limiting modes, as in the rigid- waU setting [39j. 

Definition 4.3 (Neutral limiting modes). A triple ((jjs, cts, Cg) with as positive and 
Cs real is called a neutral limiting mode if it is the limit of a sequence of unstable 
solutions {(0fc, a/c, Cfc)}fc^i of /[4.1^ ~ (4^ as k ~+ oo. The convergence of 4>k to (j)s 
is in the almost- everywhere sense. For a neutral limiting mode, as is called the 
neutral limiting wave number and Cs is called the neutral limiting wave speed. 



Lemma 14.81 below establish that neutral limiting wave numbers form the bound- 
ary points of the interval of unstable wave numbers, and thus the stability investiga- 
tion of a shear flow is reduced to find all neutral limiting modes and then study the 
stability properties near them. In general, it is difficult to locate all neutral limiting 
modes. For flows in class /C"*", nonetheless, neutral limiting modes are characterized 
by the inflection value. 

Proposition 4.4. For U G /C+ a neutral limiting mode {(j)s,as,Cs) must solve 
g3p-ra with Cs = Us. 



For the proof of Proposition 14.41 we need several properties of unstable solu- 
tions. First, Howard's semicircle theorem holds true in the free-surface setting f62l 
Theorem 1]. That is, any unstable eigenvalue c = Cr + ici (q > 0) of the Rayleigh 
equation (|4.ip - (|4.2p must lie in the semicircle 

(4.9) {Cr - i(£/mi„ + f/max))' + ^ (i((7,nin - t^max))', 

where U^in = niin[o,ft] U{y) and i/max = max[o,/,.] U{y). 

The identities below are important for future considerations. 

Lemma 4.5. If (p is a solution j[ )- |7^ with c — Cr + iCi and Ci ^ then for 
any q real the identities 

(4.10) 

+ «'H' + ^^^}^^l'/'p)dy- (Reg.(c) + ^Img.(c)) |0(Mp, 



(4.11 







[\cbf + a'lcj^l^ + ^'J^^p'V l^) dy ^ (Reg,(c) - ^Im5s(c)) |0' ih)f 



hold true, where gr is defined in J^. ?[ ) and 
(4.12) gs(c) - 



g + U'{h){U{h)--cy 
Proof. We rewrite the Rayleigh equation (|4.ip as 

U" 

(4.13) -(j)" + a^(j)+- = 0. 

U — c 

Multiplication of above by (jf and integration by parts with the boundary condition 
(1121) yield 
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Its real and imaginary parts read as 

(4.14) /" [ |<^f + \4>\^ + ^3'|2'^ I'^l' ) = ^^5.(c) 10 {hf 

Jo \ \U-c\ 



and 



h 



(4.15) c,J^ ___\^\^dy = lmgricMih)\\ 

respectively. Combining (|4.14[) and (|4.15p then establishes (|4.10[) . 
Similarly, combining the real and the imaginary parts of 

TT" \ 

2 . „2| /|2 I "-^ I i|2 1 J , *f„\\^'fu\\2 



(4.16) ^\^'\^+a^\^\^ + ^^\<l>\'jdy^g:{c)\cf,'{h)\ 

leads to ((4lT|) . □ 

Our next preliminary result is an a priori ^-estimate of unstable solutions near 
a neutral limiting mode. 

Lemma 4.6. For U G /C^, let {{4'k, ct/c, Cfe)}fe°^i &e sequence of unstable solutions 
to ^.l^ - (4^ such that ||0/c||l2 = 1- // a/c — > > and Imcfe ^ 0+ as fc —> oo 
</ien ||0fe||^2 ^ C, where C > is independent of k. 

Proof. By the semicircle theorem (j4.9p . e [C/min, t^max] for any k > 1. Thus 
Ck Cs [Umin, C^max], as A: — *■ cxD. The proof is divided into to the following two 
cases. 

Case 1: U{h) ^ Cs- The proof is similar to that of 39, Lemma 3.7]. It is 
straightforward to see that 

Imgrick) 



(4.17) |Re<7.(cfe) 



=^ Co, 



Imck 

where Cq > is independent of k. For simplicity of notations, the subscript k is 
suppressed in the estimates below and C is used to denote generic constants which 
are independent of fc. Let c = Cr + ici. We write (|4.10[) as 

(|0f + + K{y) d2/= (Re<?.(c) + ^Im5.(c)) 

Note that gr{cs) is well defined. Setting q = Us — 2 {Us — c^) in the above identity 
leads to 



" {U - Crf - {Us - Crf 

U-Cr\^+C^ 



K{y) ' „ ,r,\2^ \^?dy + c\m? 



h 



^ I K{ym^dy + C{eU'[,l. + -ml.). 
Jo = 

One chooses e sufficiently small to conclude that ||0||^ri ^ C. 

Next is the iJ^-estimate. Multiplication of (|4.13|) by —{(j)*)" and integration over 
the interval [0, h] yield 

(4.18) /'(|0'f + a' |0f )dy - «V*(c)l0(MP = [\r)"j^^dy. 

./n ./n U - C 
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In view of the Rayleigh equation for (j)* , the right side is written as 



, U" 

U ~c 



>dy 



h jj„ 



U-c 



U" 



U" 

U-c 



>dy 



iU") 



\U~c[ 



The real part of (|4.18p then reads as 



h Tjii 



U"{U ~ Cr) 
h 



Wdy- 



h ^jj„y 



{U"f 

|C/-C|2 

{U"f 
\U-c\ 



My 



M?dy,. 



where in the last equality we use (|4.14p . So 

(4.19) + 2c? \<i>'f + a\^\^)dy = 2c?R^gr{c)\m? 

Jo 

By with q ^ Us, 



\U-c[ 



Ml^dy 



-U"{U-Us) 
|C/-c|2 ' 



My (since - U"{U -Us) ^ 0) 



= i|i^ll 



(|0f + a^\ct^\^)dy - Regr{c) + (c, - Us) 



Inigr(c) 



where we use the bound (|iT7)l and that \\4>{h)\ ^ Co||(/)||_ffi. Therefore, by (|iTgi) 
and the ||(/>||j:^i bound, we get ||0||j:^2 < C as desired. 

Case 2: U{h) — Cs- The proof is similar to that of Case 1 except that one uses 
(f){h) — gs{c)(f)'{h) in place of <f)' {h) = gr{c)(j){h). It can be check that 



(4.20) 



\^egs{ck)\ 



Imgsick) 



Imcfe 



^Codick,U{h)), 



where Co > is independent of k and 



d(cfe, U{h)) = \Reck - U{h)\ + (Imcfe)' . 

Since U{h) = c^, it follows that d{ck, U{h)) — > as A: — > oo. The same computations 
as in Case 1 establish that 



+ alm^)dy ^ I K{y)\M^dy + Cod{ck,U{h))\cl,i{h)\' 
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and 

(4.21) rmf + 2al\cl^',f + at\cbk\')dy 

Jo 

= 2alReg.{ck)\Mh)\' + £ -^^\^k\'dy 
^ C{ I m' + al\cl,kndy + dick,Uih))\c^iih)f 




where C > is independent of k. Consequently, by choosing s small, 

UkWl^ ^ CiiUkWh + d{ck, U{h)) {hf ) < C2(\\M\l^ + d(ck, U(h))Uk\\m). 

where Ci, C2 > are independent of k. Since d{ck, U{h)) as fc — > 00 it follows 
from above that ||0fe||j:/2 ^ C. □ 

For U e /C+, if c is in the range of U then U{y) — c holds at a finite number of 
points [39l Remark 3.2], which are denoted by j/i < 2/2 < • • • < ym^- Let yo — 
and yrria+i — h. We state our last preliminary result. 

Lemma 4.7 (j39|. Lemma 3.5). Let (p satisfy {4-1^ with a positive and c in the 
range of U and let U{y) = c for y G {yi, 2/2, ■ . ■ , Vm^}- If't' 'is sectionally continuous 
on the open intervals {yj,yj+i), j = 0,1,..., mc, then (j) cannot vanish at both 
endpoints of any intervals [yj,yj^i) unless it vanishes identically on that interval. 

Proof of Proposition \4.4\ Let ((/)s, as, Cs) be a neutral limiting mode with > 
and Cs e [t^nin, f^max], and let {{4>k,cek,Ck)}kLi be a sequence of unstable solu- 
tions of (|4.ip - (|4.2p such that {(f>k,ak,Ck) converges to (0^, as, Cs) as fc ^ 00. We 
normalize the sequence by setting = 1. 

First, the result of Lemma l4.6l savs that ||'/>fc||//2 C, where C > is independent 
of k. Consequently, (f>k converges to 0s weakly in and strongly in H^. Then 
II 0s II //2 ^ C and II 0s II L2 = 1. Let yi,y2, - ■ ■ ,ym, be the roots of U{y) ~ Cg and let 
5*0 be the complement of the set of points {yi,y2, ■ ■ ■ ij/ms} in the interval [0, h]. 
Since 0^ converges to 0s uniformly in on any compact subset of So, it follows 
that 0" exists on Sq. Since (U — Ck)"^ , 4>k and their derivatives up to second order 
are uniformly bounded on any compact subset of 6*0, it follows that 0s satisfies 

(4.22) 0;' - a20s - TT^0s = 

U - Cs 

almost everywhere on [0, ft.]. Moreover, 
in case U{h) ^ Cg, and 

(4.24) 0s(ft) = O and 0s(O) 

in case U{h) = Cs. 
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Next, we claim that Cs is the inflection value Us- By Definition 14.11 U"{yj) = 
—K{yj){cs — Us) has the same sign for j — 1, - ■ ■ , ms, say positive. Let 

Es = UZ\{y&[0,h] : \y-y^\<d}. 

Clearly, Eg C 5*0 and U"{y) > for y e E's when S > sufficiently small. The 
proof is again divided into two cases. 

Case 1: U{h) ^ Cs- Since (j>s is not identically zero, Lemma [4.71 asserts that 
4's{yj) 7^ for some yj. If Cs were not an infiection value then near such a yj it 
must hold that 

(4.25) / — \<j)s\dy^l — —\<j)s\ dy 

Since 



IJ,' |2 , 2|i |2 , - f^min + 1) I , ,2 \ , 

l^fcl +ak\cl)k\ + \U -Ck^"^ ' ^ 



^ / 177 12 (l)krdy-snp— —. 



by Fatou's Lemma and (|4.25p it follows that 

liminf /' + al\^k\' + ^"^T^^Z^^^ M' ] dy ^ ^. 

On the other hand, (|4.10p with q ~ C/min ^ 1 yields 



"U'|2^ 21. 12 , U"{U-Umin + l) ,, ,2^ , 

ml +oikm\ + \u -ck\^ — ) ^ 

Imc/c / 

where Ci > is independent of k. A contradiction proves the claim. 

Case 2: U{h) — Cs- The proof is identical to that of Case 1 except that we use 
(|4.1ip in place of (|4.10p and hence is omitted. □ 

The following lemma allows us to continuate unstable modes until a neutral 
limiting wave number is reached, analogous to 39, Theorem 3.9]. 

Lemma 4.8. For U G /C+ , the set of unstable wave numbers is open, whose any 
boundary point a satisfies that —o? is an eigenvalue of the operator —j^ — K{y) 
on y & (0, h) with the boundary conditions |^.5[ )-| |^~5[ ). 

Proof. An unstable solution (0o, ao, co) of (|4.ip - (|4.2p is held fixed, with > and 
Imco > 0. For ri, r2 > 0, let us define 

In = {a > : |a — ao\ < ri} and — {c ^ : |c — co| < r2}, 

where denotes the set of complex numbers with positive imaginary part. Our 
goal is to show that for each a G Ir^, there exists c{a) £ for some r2 > and 
(pa G such that (</)„, a, c(q:)) satisfies (|4T|) - (|42)) . 

For a G (0, oo) and c G C^, let (/>i (z; a, c) and (/)2 (x; a, c) to be the solutions of 
the differential equation 

U" 

(I)" - a'^cf)- (f)^0 for ye(0,h) 

U — c 
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normalized at h, that is 



Mh) = 1, 



It is standard that c/fi and (/)2 are analytic as functions of c in . Let us consider 
a function on (0,oo) x C+, defined as 



where gr is defined in (|4.7p . Clearly, $ is analytic in c and continuous in a. Note 
that $(a,c) = if and only if the system (|4.ip - (l4.2p has an unstable solution 



Since <&(ao, cq) ~ and the zeros of an analytic function are isolated, $(ao, c) ^ 
on {|c — Col = r2}, for r2 > sufficiently small. Then, by the continuity of <i>(a, c) 
in a, we have $(a,c) ^ on {|c— co| = ^2}, when a € 1^ and ri is sufficiently 
small. Let us consider the function 



where a € 7^. Observe that N{a) counts the number of zeros of $(a,c) in Br2- 
Since N{aQ) > and N{a) is continuous as a function of a in it follows that 
N{a) > for any a G . This proves that the set of unstable wave numbers is 
open. 

By definition, the boundary points of the set of unstable wave numbers must be 
neutral limiting wave numbers, say as- Proposition 14.41 asserts that — must be 
a negative eigenvalue of the operator — K{y) on (0, h) with (|4.5P " (|4.6p . This 
completes the proof. □ 

4.2. Proof of Theorem 14. 2i The proof of Theorem 14.21 is to examine the bi- 
furcation of unstable modes from each neutral limiting mode. This reduces to 
the study of the bifurcation of zeros of the algebraic equation $(a,c) = from 
{as, Cs) as is pointed out in the proof of Lemma 14.81 However, the differentiability 
of $(q!, c) in c at a neutral limiting mode (a^jCs) can only be established on half 
of the c— neighborhood, and thus the standard implicit function theorem does not 
apply. To overcome this difficulty, we construct a contraction mapping in such a 
half neighborhood and prove the existence of an unstable mode near (as,Cs) as is 
done in |39j for the rigid-wall case. Moreover, the existence of an unstable mode 
can only be established for wave numbers slightly to the left of as- Therefore, 
unstable modes with wave numbers slightly to the left of Ug can be continuated to 
the zero wave number. Our method is almost the same as in the rigid-wall setting 
([31 Section 4]). 

Below we construct unstable solutions for wave numbers slightly to the left of a 
neutral limiting wave number. 

Proposition 4.9. For U (y) in the class IC^ , let Us be the inflection value and 
2/1: 2/2, ■ • ■ , J/nis be the inflection points, as such U{yj) — Us for i = 1,2,..., to^. // 

{(ps, as,Us) is a neutral limiting mode with as positive, then for e G (eo,0), where 
|eo| is sufficiently small, there exist and c(e) such that 



(4.26) 



$(a,c) = (f>i{0-a,c) + gr{c)(l)2{0;a,c), 



4>a{y) = <j)i{y;a,c) + gr{c)(j)2{y;a,c). 




d^/dc{a, c) 
<i>(a, c) 



(4.27) 
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(4.28) (j)',ih) ^ griUs + c{e))^,ih) and 0^(0) = 
with Inic(£) > 0. Moreover, 

(4.29) lim c(e) = 0, 

£^0 — 

(4.30) 



here p.v. means the Cauchy principal part and 

2g j_ U'{h) 

(4.31) A 



(U(h)-U,Y + (U{h)-Us)^ U{h) ^ Us 

if U{h) = C/, 



Proof. As in the proof of Lemma [4.81 for c e C+ and e < 0, let (j)i{y;e,c) and 
02 (y; £, c) be the solutions of 

(4.32) ^" -(al+e)cl>~-—^ = for y e {0,h) 

U - Us - c 

normalized at h, that is 

Ul{h) = l, 02(M=O, 

\0'i(M = O, <P'^{h) = l. 

It is standard that 0i and 02 are analytic as a function of c in and that 0i 
and 02 are linearly independent with Wronskian 1. The neutral limiting mode is 
normalized so that 0s (ft.) = 1 and 0s (ft) — gr{Us). The proof is again divided into 
two cases. 

Case 1: U{h) ^ Us- Let us define 

0o(y;£,c) = (j>i{y;e,c) + gr{Us + c)(l)2{y;e,c), 
$(e, c) = 01 (0; £, c) + g^(C/s + c)02(O; e, c), 

where gr is given in (|4.7p . It is readily seen that 0o solves (I4.32p with 0o(ft) = 1 
and 0o(ft) = gr{Us + c). It is easy to sec that $(e, c) is analytic in c S C+ and 
differentiable in e. Note that an unstable solution to (|4.27p - (|4.28p exists if and only 
if <&(£,£) = for some Imc > 0. The Green's function of (I4.32p is written as 

G{y,y';e,c) = 0i(y; e, c)02(2/'; e, c) - 02(2/; e, c)0i(y'; e, c), 

where 0i(2/;e, c) is the solution of (|4.32p with 0i(ft) = 1 and 0'i(ft) = gr{Us)- A 
similar computation as in pp. 336] for (l)j{y;e,c) (j — 1,2) yields that 

(4.33) 'de^'^'^^^^J 

and 
(4.34) 

9$ -U" d 

— (e,c) = j G(y,0;e,c)^^^— — y2 0o(y;£,c)dy + — 5,.(?7s + c)02(O;e,c). 

Let us define the triangle in C+ as 

^{R,b) = {c-r + iCi ■■ \Cr\ < RCi , < Q < 5} 
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and the Cartesian product in (0,oo) x C+ as 

E{R,bi,b2) = (-^2,0) X A(fl_f,j), 

where i?, 61, 62 > are to be determined later. 
We claim that: 

(a) For fixed R, both (pi (•; e, c) and (^o('; c) uniformly converge to (j)^ in [0, h] 
as (£,c) (0,0) in -^(fl^bj.bj). That is, for any S > there exists some &o > such 
that whenever 61,62 < 60 and (e, c) € £'(fl,6i.b2) the inequalities 

pi(-;e,c) - 0s||ci ' \\^oi-;£,c) - =^ 

hold. 

(b) (j)2{-]e,c) converges uniformly to (j)2 {y;0,0) in the sense of (a). We denote 
02 (y;0,0) = (j)z (y) , then 0^ (0) = —-^pjoj since the Wronskian of (^i,(/'2) and its 
limit (0s, (/iz) is 1- The proof of (a) and (b) is very similar to [321 PP- 337] and we 
skip it. 

In the appendix, we prove that 

uniformly as e — > 0— and c — > in i?(fl,bi,f,2) ; where A is defined by (|4.3ip and yj 
(j = 1, . . . , TOs) are the inflection points for Us- Let us denote 

1 



Lemma |4 . 71 asserts that 0^ is nonzero at one of the inflection points, say (j)s{yj) ^ 0. 
Note that by [SS Remark 4.2], for class K.+ flows, U'{y.j) 7^ for j = 1, 2, ... , m^. 
Consequently, D <Q. 

The remainder of the proof is identically the same as that of j39» Theorem 4.1] 
and hence we only sketch it. Define 

(4.38) /(£,c) = $(e,c)-Be-(C + r'i)c, 

(4.39) ^(^'^)— 

G + lU G + iD 

Note that for each e < fixed, a zero of $(e, •) corresponds to a fixed point of the 
mapping F{e,-). It is shown in [39, pp. 338-339] that for e G (eo,0), where Jeoj is 
sufficiently small, the mapping i^(e,-) is contracting on A(j^ b(j)) for some i? > 
and 

6(e) = ~2DB{C^ + D^y^e. 
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So for each e G (eo,0) there exists a unique c{e) £ A(fl^fc(£)) such that F{e,c{e)) = 
c{e) and thus <i>(£,c(£)) = 0. It can be also shown that c(e) is diffcrentiable in 
e in the interval (eo,0). Since c(e) G A(fl it is immediate that (I4.29P holds. 
Finally, differentiation of <i>(e,c(e)) = yields 

which, in view of ((05)) and (HSU), imphes ((Oil)) . 

Case 2: ?7(/i) = C/g. The proof is almost identical to that of Case 1, and thus we 
only indicate some differences. Let us define 

$(£, c) = gs{Us + c)(f>i (0; e, c) + 02 (0; e, c), 

where gs is defined in (|4.12|) . Note that (j)s{h) — 0s (0) = 0. Thereby, the Green's 
function is written as 

G{y,y';e,c) = 4>i{y;£,c)(f>2{y' ■,€,€) - 02(2/; e, c)0i(2/'; e, c). 

The same computations as in Case 1 yield that 

a$ 1 



de 0',(O) 



tdy 



and 



uniformly as £ — s- 0— and c — s- in E(^]i i,^ i,^y This completes the proof. □ 



Proof of Theorem \4.2\ Let —ajf < —c^n-i < ■ ■ ■ < —a{ < be negative eigenval- 
ues of the operator — ^ — K{y) on {Q,h) with boundary conditions (I4.5|) - (I4.6|) . 
That is, a AT — Ofmax, where — a^j^x is defined either in (I4.8ap or ()4.8bp . We de- 
duce from Lemma [4.81 and Proposition 14.91 that to each a e (0, aA?) with a ^ aj 
(j = 1, . . . , N) an unstable solution is associated. Our goal is to show the instability 
at a = aj for each j = 1, . . . , iV — 1. 

Case 1: U{h) ^ Us- Let {(0fe, afe, Cfe)}fe°;i be a sequence of unstable solutions 
such that afe — > Oij+ as /c — > 00. After normalization, we may assume (t>kih) = 1. 
Note that </>/. satisfies 

(4.40) ^'^-al(f>k-j^^<f>k^O for y G (0, /i) 

and <t>'k{h) = gr{ck), 4>k{Q) = 0. Below we will prove that Imcfc ^ (5 > 0, where 5 is 
independent of fc. Since the coefficients of (|4.40p and gr{ck) are bounded uniformly 
for k, the solutions 4>k of the above Rayleigh equations are uniformly bounded in 
C^, and subsequently, (j)k converges in as A: ^ 00, say to 0oo- The semicircle 
theorem (|4.9p ensures that Cfe ^ Coo as fc — *■ 00. Note that Imcoo ^ ^ > 0. By 
continuity, (j)oo satisfies 

U" 

'i^L - "j'/'oo - TT^ 0OO = for ?/ e (0, /l), 

with 0oo(ft-) = 1, 0Jx3(^) = .9r(coo) and 0oo(O) = 0. That is, (0oo,aj,Coo) is an 
unstable solution of (lOl-dOl. 
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It remains to show that {Imcfc} has a positive lower bound. Suppose on the 
contrary that Imcfc — > as fc —> cx3. 

We claim that ||0fc||i2 ^ C, where C > is independent of k. Otherwise, 
||(/)fc||i2 ^ oo as fc ^ oo. Let ipk = (j^k/WML^, then ||'^fc||L2 = 1. Lemma [Jj] 
then dictates that ||(y9fc||^2 ^ C independently of k. Subsequently, Proposition 
14.41 ensures that {(pk,cxk,Ck) converges to a neutral limiting mode {cps, as,Us)- By 
continuity, ||<^s||l2 = 1 and 

^'^ - a]ifs + K{y)Lps =0 for y e (0, h). 

On the other hand, tfs{h) = f'si^) = since (pk{h) = I/II'/'/cIIl^ and ifikW ~ 
5r(cfc)/||0A:||L2 — > as fc — !■ oo. Correspondingly, ips = on [0, ft.]. A contradiction 
proves the claim. 

Since ||0fc||L2 is bounded uniformly for k, Lemma 14.61 and Proposition 14.41 applv 
and ||0fc||_ff2 ^ C, Ck Us and (f>k — > 0s in C^, where (ps satisfies 

- - TT^^s = for y e (0, h) 

U - Us 

with 4)s{h) — 1, 4)'s{h) — gr{Us) and (t>s{^) ~ 0. An integration by parts yields that 

= ^ l^^sicl^l - aUk - YJ^^M - M€ - c^jc^s - ^7^-^^)) dy 

= {aj-al) (jJsh dy - (ck - Us) — -— — —r(t)s(l)kdy + grick) - gr{Us). 

Jo Jo {U~ck){U~Us) 

Let us denote 

Bk^ 4>s(l>kdy, 



n / ^" A. A A . 9r{ck) - gr{Us) 

B>k = - tTt T777 TTT9s<Pkdy + 



/o {U ~Ck){U -UsY""" " Ck~Us 
It is immediate that limfc^oo Bk — \(t>s\'^dy. We shall show in the appendix that 



(4.41) lim I?,=A + ^7^g^M02( 

fc^oo ^ \U{yj)\ 



where A is defined by (|4.3ip and Uj (j = 1, 2, . . . , nis) are inflection points corre- 
sponding to the inflection value Us- Since Im (limfe_»oo -D^) > (see the proof of 
Proposition l4.9p it follows that 

Imcfe = {al - a^)Im{Bk/Dk) < 

for k large. A contradiction proves that {Imc^} has a positive lower bound, uni- 
formly for k. 

Case 2: U{h) = Us- We normalize 4>k so that 0jj(ft.) = 1 and (pkih) — gs{ck)- 
The proof is identically the same as that of Case 1 except that 

U" , ^ ^ , gsick) 

Dk = - TTt T777 TTT^sA dy -' 



/o {u~ck)iu-Usy'^^ " ck-Us' 

We shall show in the appendix that 

(4.42) ^hm^i.. ^ -| ^^^(^.) f u^fs<^y- 
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This proves that there exists an unstable solution for each a G (0, ctmax)- 

It remains to prove linear stability in case either the operator — K{y) on 

y G (0, h) with ()4.5|) - ()4.6|) is nonnegative or a ^ amax- Suppose otherwise, there 
exists an unstable mode at a wave number a ^ amax- By Lemma 14.81 can 
continuate this unstable mode for wave numbers larger than a until the growth 
rate becomes zero. By Lemma |6.3[ this continuation must stop at a wave number 
as > a, where there is a neutral limiting mode. Then by Proposition 14. 4[ —a1 is a 



negative eigenvalue of — K (y) with (|4.5p - (|4.6p . But —a1 < — an^ax which is a 
contradiction to the fact that — a^^^x the lowest eigenvalue of — K (y) with 
This completes the proof. □ 

Remark 4.10. For U G /C+, let —a^ be the lowest eigenvalue of — ^ — K{'y) on 
y £ (0, /i) with the Dirichlet boundary conditions 4>{h) = = (/)(0). If U{h) = Us 
then ttmax = Oid- We claim that if U(h) ^ Us then a,„ax > otd- To see this, let 
be the eigenfunction of — ^ — K{y) on y £ (0, h) with the Dirichlet boundary 
conditions corresponding to — and let be the eigenfunction of —-j^ ~ K{y) 
with the boundary conditions (|4.5p - (|4.6p corresponding to — a^^x- By Sturm's 
theory, we can assume 4>di 4>m > on y € (0, /i). An integration by parts yields that 

= / - "max'/'™ + K{y)<j),n) " 0m (0d " "d^d + K{y)(j)d)) dy 

Jo 

= -(l)m,{h)4''d{h) + {al - "iax) / Mmdy. 



JO 

Since (f>'d{h) < the claim follows. 

For flows in class /C+, the lowest eigenvalue of (|4.4p measures the range of in- 
stability, for both the free surface and the rigid wall cases. That means, (0, amax) 
is the interval of unstable wave numbers in the free-surface setting (Theorem 14. 2p 
and (0, ad) is the interval of unstable wave numbers in the rigid-wall setting j39l 
Theorem 1.2]. The fact that amax > ctd thus indicates that the free surface has a 
destabilizing effect. 

4.3. Monotone unstable shear flows. In general, for a given shear flow profile 
in the class /C+, one might show the existence of a neutral limiting mode and thus 
the existence of growing modes, by numerically computing the negativity of (|4.8p . 
For monotone flows with one inflection point, however, the existence of neutral 
limiting modes can be shown ([52]) by a comparison argument. 

Lemma 4.11. For any monotone shear flow U (y) with exactly one inflection point 
ys in the interior, there exists a neutral limiting mode. That is, has a 

nontrivial solution for which c = U{ys) and a > 0. 

Proof. This result is given in Theorem 4 of 62J. Here we present a detailed proof 
for completeness and also for clarification of some arguments in |62| . 

We consider an increasing fiow U (y) only. A decreasing fiow can be treated in 
the same way. Let Us = U{ys) be the infiection value. Denoted by 0q the solution 
of the Rayleigh equation 

Yr~U - = for ye{0,h) 
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with (j)a{Q) ~ and (p'^iO) — 1- As in the proof of Lemma [231 an integration of the 
above over (0, h) yields that 

{U{h) ~ Us)^'^{h) ~ (C/(0) - Us) ~ Mh)U'{h) - ([/ - Us)c^ady - 0, 

Jo 

and thus 

U{0)-Us U'jh) [\u^JJU d 

Mh) iU{h) - Us)Mh) U{h) ^Us {U{h) - Us)Mh) Jo 
It is straightforward to see that the boundary condition (|4.2p is satisfied if and only 
if the function 

/(") - nr^^^ rr ^1 + TTTTT^ FTTTTTTT / - Us)^o.dy ■ 



{u{h) - Us)Mh) iu{h) ~ Us)Mh) Jo {Ue - u{h)f ■ 

vanishes at some a > 0. 

We claim that (paiu) > on y G (0, ft-] for any a ^ 0. Suppose otherwise, let 
ya G (0, h] to be the first zero of 0q other than 0, that is, (j)a{ya) = and (paiy) > 
for y G (0,?/q,). Then, y^ > ys must hold. Indeed, if =^ 2/s were to be true, then 
Sturm's first comparison theorem would apply to (j)a and U — Us on [0, ya] to assert 
that U — Us must vanish somewhere in (0, ya) C (0, ys). This uses that 

{U-Usf + j^{U-Us)^i). 

A contradiction then asserts that ya > ys- Correspondingly, <f>a and U — Us have 
exactly one zero in [0,ya]. On the other hand, by Sturm's second comparison 
theorem [3T], it follows that 

<(ya) > U'{ya) 



(t)a{ya) U{ya)-Us' 

This contradicts since (j)a(jja) — and the left hand side is —oo. Therefore, <f>a{y) > 
for y E {0,h] and for any a ^ 0. In particular, (paih) > for any a ^ 0. 
Consequently, / is a continuous function of a and /(O) < 0. 

It remains to show that /(a) > for a > big enough. Thereby, by continuity 

^ M and > M. Let us denote by 4>i and 



u" 



Us-U 



/ vanishes at some a > 0. Let 
4>2 the solutions of 

01 + (M- a^)^! = and (j)^ + {-M - a^)4>2 = for ye(0,/i), 
respectively, with (j)^ (0) — and (j)^ (0) = 1 . It is straightforward that 
1 . , 7T . , , s 1 



^i(y) = / sinh - My and (?!)2(y) = . „ sinh Va^ + My. 

Va^ — M Va^ + M 

As in the proof of Lemma [2. 71 (a). Sturm's second comparison theorem [31] implies 
that 

: sinh \/ cP- — M y ^ <i>a {y) ^ — sinh \/ a'^ + My. 



\Ja^ - M ' Vo^ + M 

This together with the monotone property of U establishes that /(a) ^ Cia — C2 
for some constants Ci,C2 > 0. For details we refer to [62l Theorem 4]. Thus, 
/(a) > if Of > is sufficiently large. This completes the proof . □ 

Since a monotone flow with one inflection value is in class JC^ , the above lemma 
combined with Theorem 14. 21 asserts its instability. 
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Corollary 4.12. Any monotone shear flow with exactly one inflection point in 
the interior is unstable in the free-surface setting, for wave numbers in an interval 

(0, amax) with amax > 0. 

Remark 4.13. In the free-surface setting, there are two different kinds of neutral 
modes, which are solutions to the Rayleigh system ()4.1|) - (I4.2|) with Imc — 0. Neu- 
tral limiting modes have their phase speed in the range of the shear profile. For 
flows in class /C+, moreover, the phase speed of a neutral limiting mode must be 
the inflection value of the shear profile (Proposition [4]4]) and it is contiguous to un- 
stable modes fProposition I4.9[) . On the other hand. Lemma [2.51 shows that under 
the conditions U"{h) < and U{h) > U{y) for h y, a, neutral mode exists with 
the phase speed c > max U. Such a neutral mode is used in the local bifurcation 
of nontrivial periodic waves in Theorems 12.21 and 12.61 In view of the semicircle 
theorem (j4.9p . however, such a neutral mode is not contiguous to unstable modes. 
This implies that neutral modes governing the stability property are different from 
those governing the bifurcation of nontrivial waves. This is an important difference 
in the free-surface setting. Since in the rigid wall setting, for any possible neutral 
modes the phase speed must lie in the range of U, which follows easily from Sturm's 
first comparison theorem. For class /C^ flows, such neutral modes are contiguous 
to unstable modes ([39]) and the bifurcation of nontrivial waves from these neutral 
modes can also be shown ([9]). Thus, in the rigid- wall setting, the same neutral 
modes govern both stability and bifurcation. 

Remark 4.14. In [TH], the J'-formal stability was introduced via a quadratic form, 
which is related to the local bifurcation of nontrivial waves in the transformed 
variables (see also Introduction), and it was concluded that this formal stability of 
the trivial solutions switches exactly at the bifurcation point. Below we discuss two 
examples for which the linear stability property does not change along the line of 
trivial solutions passing the bifurcation point, which indicates that the j7-formal 
stability in [19] is unrelated to linear stability of the physical water wave problem. 
However, the j7-formal stability results [19] do give more information about the 
structure of the periodic water wave branch. 

Let us consider a monotone increasing flow U {y) ony G [0, h] with one inflection 
point ys G (0,/i), for example, U{y) = as'mb{y — h/2) on y G [0,h] for which 
ys = h/2. By lemma H. Ill and Corollarv l4.121 such a shear flow is equipped with a 
neutral limiting mode with c ^ U (jjs) and a = amax > 0, and it is linearly unstable 
for any wave number a € (0, amax)- In addition, by lemma [2.51 and Theorem 
12.61 for any wave number a € (0, ckmax) this shear flow has a neutral mode with 
c(a) > U{h) = maxJJ. Moreover, such a neutral mode is a nontrivial solution to the 
bifurcation equation (|2.9p - (|2.10p . and thus there exists a local curve of bifurcation 
of nontrivial waves with a wave speed c{a) and period 2T:/a. Let po and 7 be the 
flux and vorticity relation determined by U (y) , c(a) and h via (|2.1ip . Consider the 
trivial solutions with shear flows U{y]fj,) defined in Lemma |2. 11 with above po, 7 
and the parameter ^. The bifurcation point U{y; /io) = U [y) — c{a) corresponds 
to jio = {U{h) — c{a))'^ . The instability of U{y;^o) at the wave number a is 
continuated to shear flows U (y; /i) with fi near /xq, which can be shown by a similar 
argument as in the proof of Lemma [4?8l So at the bifurcation point ^0, there is NO 
switch of stability of trivial solutions. 

Let us consider U € C^([0, /i]) satisfying that U'{y) > and U"{y) < in 
y & [0,h]. For such a shear flow Lemma [2.51 and Theorem 12.61 applies as well and 
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there exists a local curve of bifurcation of nontrivial waves for any wave number 
a which travel at the speed c(a) > max U, where c{a) is chosen so that the bifur- 
cation equation (|2.9p - (|2.10p is solvable. For this bifurcation flow U {y) — c{a), the 
vorticity relation 7 determined via (|2.1ip is monotone since U" does not change 
sign. Consequently, any shear flows U{y;fj,) defined in Lemma 1 2 . II with the same 
7 has no inflection points, as also remarked at the end of Section 3. Therefore, 
by Theorem 16.41 all trivial solutions corresponding to these shear flows U(jj;ii) are 
stable. This again shows that the bifurcation of nontrivial periodic waves does not 
involve the switch of stability of trivial solutions. 

5. Linear instability of periodic water waves with free surface 

We now turn to investigating the linear instability of periodic traveling waves 
near an unstable background shear flow. Suppose that a shear flow ([/(j/),0) with 
U E C'^^^ {[0, ho]) and U E /C+ has an unstable wave number a > 0, that is, 
for such a wave number a > the Rayleigh system (|4.1|) - (|4.2|) has a nontrivial 
solution with Imc > 0. Suppose moreover that for the unstable wave number a > 
the bifurcation equation l|2.9p - (|2.10p is solvable with some c(a) > maxU. Then 
Remark l2.4l and Theorem 12 . 2 1 applv to state that there exists a one-parameter curve 
of small-amplitude traveling- wave solutions (j]e(x),'ipeix,y)) satisfying (|2.2p with 
the period 27r/a and the wave speed c(a), where e ^ is the amplitude parameter. 
A natural question is: are these small-amplitude nontrivial periodic waves generated 
over the unstable shear flow also unstable? The answer is YES under some technical 
assumptions, which is the subject of the forthcoming investigation. 

5.1. The main theorem and examples. We prove the linear instability of the 
steady periodic water-waves (?7e(a;), ipeix, y)) by finding a growing-mode solution to 
the linearized water-wave problem. As in Section [H let 

Pe = {ix,y) : < X < 2TT/a, <y < 7]f_{x)} and = {{x,ri^{x)) : < x < 27r/a} 

denote, respectively, the fluid domain of the steady wave {rj^{x),il)^{x,y)) of one 
period and the steady surface. The growing-mode problem (|3.9p of the linearized 
periodic water-wave problem around (ri^{x),'ijje{x,y)) reduces to 

(5.1a) A7A + 7'(V',)V'-7'(^e) / Xe^'7/j{X,{s),Y,{s))ds = in P,; 



(5.1b) \r]{x) + -^(■0ey(x,77,(x))?7(a;)) = --^-^(a;, 77^(2;)); 

(5.1c) P{x, f],{x)) + Pev{x)v{x) = 0; 

(5. Id) X-tpnix) + -^{■ip^y{x,ri^{x))'>pnix)) = --^P{x,'q^{x)) - Vt-^^{x,r]e{x)); 
(5.1e) V(a;,0) = 0. 

Here and in sequel, let us abuse notation and denote that P^y{x) = P^y{x,ri^{x)), 
that is, the restriction of P^y{x,y) on the steady wave-profile y — tjeix). By Theo- 
rem [221 Pey{x) = Pey{x,7je{x)) = -g + 0(e). Rccall that 

ipnix) = dyil^{x,rie{x)) - ■qex{x)dxil!{x,rie{x)) 

is the derivative of il^ix, '7e(x)) in the direction normal to the free surface (x, rj^ix)) 
and that SI = 7(0) is the vorticity of the steady flow of ■ipt{x,y) on the steady 
wave-profile y — ri^{x). Note that is a constant independent of e. 
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Theorem 5.1 (Linear instability of small-amplitude periodic water-waves). Let 
the shear flow U [y) G C^+'^([0, /iq]); P G (Oj 1); in class JC^ . Suppose that 
U{h()) 7^ Us, where Us is the inflection value of U , and that amax defined by ( [^.(^a[ ) 
is positive, as such Theorem \4-.S\ applies to find the interval of unstable wave numbers 
(0, Q;,nax)- Suppose moreover that for some a € (amax/2, ckmax) there exists c(a) > 
max[/ such that the bifurcation equation 12.9\) - i2. 1 0|) has a nontrivial solution. 
Let us denote by {ri^{x),t}j^{x,y)) the family of nontrivial waves with the period 
2-K/a and the wave speed c{a), bifurcating from the trivial solution rio(x) = Hq 
and {ipQyiy), —ijjoxiy)) = U{y) — 0(0;), 0), where e ^ is the amplitude parameter. 
Provided that 

(5.2) g + U'{ho){Uiho)-Us)>0, 

then for each e > .sufficiently small, there exists an exponentially growing solu- 
tion {e'^*ri{x),e'^*'4'{x,y)) of the linearized system US. 6]) . where ReA > 0, with the 
regularity property 

{V{x),i^{x,y)) e C'+P{[0,27r/a]) x C'+f'iA). 

Remark 5.2 (Examples). As is discussed in Remark 14.141 any increasing flow shear 
flow U € C^+^([0, /lo]), /3 G (0,1), with exactly one inflection point in y e {0,ho) 
satisfies amax > 0. Moreover, Lemma 12.51 applies and small-amplitude periodic 
waves bifurcate at any wave number a > 0. Since (j5.2p holds true, therefore, by 
Theorem 15.11 small- amplitude periodic waves bifurcating from such a shear flow at 
any wave number a G (Q;niax/2, amax) are unstable. Below, we discuss in details 
such an example: 

(5.3) U{y) = asmb{y-ho/2) for ye[0,/io], 

where hQ,b > satisfy hob ^ tt and a > is arbitrary. 

(1) The shear flow in ()5.3|) is unstable under periodic perturbations of a wave 
number a G (0,amax), where amax > 0. Note that in the rigid- wall setting [39], the 
same shear flow is stable under perturbations of any wave number. This indicates 
that the free surface has a destabilizing effect. This serves as an example of Remark 
14.101 since amax > and ad = 0. 

(2) The amplitude a and the depth ho in (|5.3p may be chosen arbitrarily small, 
and the shear flow as well as the nontrivial periodic waves near the shear flow are 
unstable for any wave number a G (0, amax), which contrasts with the result in [12] 
that small-amplitude rotational periodic water-waves are jy-formally stable if the 
vorticity strength and the depth are sufficiently small. Thus, as also commented 
in Remark ()4.14p , the j7-formal stability in is not directly related to the linear 
stability of water waves. Indeed, while dj{r], tp) = gives the equations for steady 
steady waves, the linearized water-wave problem is not in the form 

which is implicitly required in |19^ in order to apply the Crandall-Rabanowitz theory 
[H] of the exchange of stability. 

(3) Our example (j5.3p also indicates that adding an arbitrarily small vorticity 
to the irrotational water wave system of an arbitrary depth may induce instability. 
That means, although small irrotational periodic waves are found to be stable 
under perturbations of the same period [48], [55], they are not structurally stable; 
Vorticity has a subtle influence on the stability of water waves. 
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The proof of Theorem 15.11 uses a perturbation argument. At e = the trivial 
sohition {rio{x),'ipo{x,y)) corresponds to the shear flow {U{y) — 0(0;), 0) under the 
flat surface {y = ho}. To simphfy notations, in the remainder of this section, we 
write U{y) for U{y) — c{a), as is done in Section[5l Thereby, U{y) < 0. Since a is an 
unstable wave number of U (y) , there exist an unstable solution (f>a to the Rayleigh 
system (j4.1|) - (j4.2|) and an unstable phase speed Cq.. That is, (pa ^ 0, IuiCq, > and 

U" 

< - + 77 0a ==0 for y e (0, ho), 

U - Ca 

-^a «o) = 77777-1 T2+ 777771 jMho), 

\[U[ho)~Cay U{ho)-CaJ 

0a (0) =0. 

This corresponds to a growing mode solution satisfying (|5.ip at e = 0, where Aq = 
—iaCa has a positive real part. Our goal is to show that for e > sufflciently small, 
there exists near Aq such that the growing-mode problem (|5.ip at {rj^(x), 'ipe(x, y)) 
is solvable. First, the system (15. ip is reduced to an operator equation deflncd in a 
function space independent of e. Then, by showing the continuity of this operator 
with respect to the small-amplitude parameter e, the continuation of the unstable 
mode follows from the eigenvalue perturbation theory of operators. 

5.2. Reduction to an operator equation. The purpose of this subsection is to 
reduce the growing mode system (15. ip to an operator equation on Lpp^(iSc). Here 
and elsewhere the subscript per denotes the periodicity in the x-variable. The idea 
is to express rj{x) on and ipi^iU) (and hence P{x,r]^{x))) in terms of 

tp{x,'n,{x)). 

Our flrst task is to relate r]{x) with ^p{x, ri^{x)). 

Lemma 5.3. For |A— AqI ^ (ReAo)/2, where Aq — —iaCa, let us define the operator 
: LlJS,) Ll,,{S,) by 

C'cl>{x) = --4t0(^) + , , ! r Xe'^^^'^i;-y\x')^{x')dx' 

e^"(-')^-i(x')0(x')dx', 



■0£j^(a;)e^'^(^) (1 - e^°(27r/a)^ 



where a{x) = Jq ^~y^ {x' , r]^{x'))dx' . For simplicity, here and in the sequel we iden- 
tify ipeyix) with ipeyix^rj^ix)) and 4>(x) with (j3{x,rj^(x)), etc. Then, 

(a) The operator is analytic in A for |A — Ao| ^ (ReAo)/2, and the estimate 

holds, where K > is independent of X and e. 

(b) For any (p G ^^^^{3^), the function ip — C^(j) is the unique Lp^^{Se)- weak 
solution of the first-order ordinary differential equation 

(5.6) X^p + -^{ip^y{x)ip) ^ --^(j). 

If, in addition, £ Cpg^(5c) then (p g Cp^^{Sf) is the unique classical solution of 

(TIT 
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Proof. Assertions of (a) follow immediately since tpf^ylx) < and thus a{x) < and 
since ReA ^ RcAo/2 > 0. 

(b) First, we consider the case G Cp^^.{S^) to motivate the definition of C'*'. 

Let us write (|5.6p as the first-order ordinary differential equation 



dx^ ip^y \ dx j ijj^y dx ' 
which has an unique 27r/Q:-periodic solution 



1 /•2Tr/a , 



An integration by parts of the above formula yields that f{x) = C^(j) is as defined 
in (f53|) . 

In case <j) G -^pcrl'^e); the integral representation ()5.5p makes sense and solves 
(|5.6p in the weak sense. Indeed, an integration by parts shows that (p defined by 
(|5.5p satisfies the weak form of equation (|5.6p 

27r/a / d d \ 

\ip{x)h{x) — tpfy{x)(p{x) — h(x) — 4>(x) — h(x) ] dx = 
dx dx J 

for any 27r/a-periodic function h g 7Jppj.([0, 27r/a]). In order to show the unique- 
ness, suppose that (p g Lp^^{Se) is another weak solution of (|5.6p . Let (pi — (p — (p. 
Then, (pi G Lp^^.{S^) is a weak solution of the homogeneous differential equation 

A<Pl + ^ {lpey{x)(pi) = 0. 

It is readily seen that Jq^^" (pi{x)dx — 0. Note that hi{x) — J^^ (pi{x')dx' defines a 
27r/Q!-periodic function in Hp^j.{[0, 2-K/a]). A/Iultiplication of the above homogeneous 
equation by (A/ii)* and an integration by parts then yield that 

27r/Q 

2 



= Re/ [\\\' ^i{x)hl{x)-^,y{x)\*^i{x)[—hi{x)] ]dx 



|A|' j^j'"^ ^ (i \h,\^) dx - RcA j^j'"^ ^,y{x) l^il' dx 

2-K/a 



= — ReA / ^ty{x)\Lpi^ dx. 
Jo 

Here and elsewhere, the asterisk denotes the complex conjugation. Since ReA > 
and ipey{x) < 0, it follows that tpi = 0, and in turn, ip = ip. This completes the 
proof. □ 

Formally, the operator can be written as 

CV(x) = -U+^ {i^,yix)^ix))\ j-c^{x). 
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Let US denote f{x) = 'ip{x,r]^{x)). With the use of then the boundary condi- 
tions (|5.1b )- (|5.1d ) are written in terms of / as 

rj{x) - C\f{x), 

Pix,f^,ix)) = -P,y{x)C^f{x), 

= -C^{Pcy{x)C^ + nid)f{x), 

where id : L'^{Se) — > L^{Se) is the identity operator. 

Our next task is to relate ip{x,y) in with f{x) — 'ijj{x,ri^(x)). Given b G 
Lp^j.{Se), let ipb G H^iPt) be a weak solution of the elliptic partial differential 
equation 

.0 

(5.7a) A^ + 7'(V.)^-7'(V'.) / \e^'i>{X,{s),Y,{s))ds ^ Q in V, 



(5.7b) ipriix) := dyip{x,r]^{x)) - ri^x{x)d^ip{x,ri^{x)) = b{x) on 5^, 

(5.7c) V(a;,0)=0 

such that ^i, is 27r/a-periodic in the a;- variable. Lemma 15.61 below proves that the 
boundary value problem (|5.7p is uniquely solvable and ipi, G H^{'D^) provided that 
|A — Aol < (ReAo)/2 and {r]^{x),ipe{x,y)) is near the trivial solution with the flat 
surface y — ho and the unstable shear flow {U{y),0) given in Theorem 15.11 This, 
together with the trace theorem, allows us to define an operator 7^ : Lp^j.{Se) 
LleriS.) by 

(5.8) %b{x) = ipbix, Veix)), 

which is the unique solution iph of (|5.7p restricted on the steady surface S^. 
Be definition, it follows that 

f{x) = lfj{x, ?7e(x)) = %tpn{x). 

This, together with the boundary conditions written in terms of / as above yields 
that 

(5.9) f^-%C\P,y{x)C^ + nid)f. 

The growing- mode problem (j5.ip is thus reduced to find a nontrivial solution f{x) — 
tp{x,T]e{x)) G Lp^^.{Se) of the equation (|5.9|) . or equivalently, to show that the 



operator 



id + %C^{P,y{x)C^ + md) 



has a nontrivial kernel for some A G C with Re A > 0. 

The remainder of this subsection concerns with the unique solvability of (j5.7p . 
Our first task is to compare (|5.7p at e = with the Rayleigh system, which is useful 
in later consideration. 

Lemma 5.4. For U [y) in class /C+, y G [0, h^] and U{ho) ^ Us, let — a^ax 
lowest eigenvalue of ~ K{y) for y G (0, Hq) subject to the boundary conditions 

where K is defined in Let — be the lowest eigenvalue of —-^ — K{y) for 

y G (0, /lo) subject to the boundary conditions 

(5.11) (t)'{ho) = and 0(0) = 0. 
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Ifg + U'{hQ){U{ho) - Us) > 0, then amax > a«. 

Proof. The argument is nearly identical to that in Remark 14.101 Let us denote by 
0m the eigenfunction of — ^ — K{y) on y e (0, /iq) with (|5.10p corresponding to 

the eigenvalue —ct^g^^ and by (/)„ the eigenfunction of — ^p- — K(i/) on j/ G (0, Hq) 
with (|5.1ip corresponding to the eigenvalue —a^- By the standard theory of Sturm- 
Liouville operators we can assume 4>m > and 0„ > on y g (0, Hq). An integration 
by parts yields that 











- "max) / (Pn4>mdy A . , . — j (/>„i (fto)'?>n (fto) • 



The assumption g + U' {ho){U{ho) — Us) > then proves the assertion. □ 
Our next task is the unique solvability of the homogeneous problem of ((57 



Lemma 5.5. Assume that C/(/io) ^ Us and g + U'{ho){U{ho) - Us) > 0. For e > 
sufficiently small and \X — Ao| ^ (ReAo)/2, the following elliptic partial differential 
equation 

M 

(5.12a) AV' + 7'(V'e)V'-7'(V'e) / \e^'ip{X,{s),Y,{s))ds ^ Q inV^ 

J — oo 

subject to 

(5.12b) ^{x,0)=0 

and the Neumann boundary condition 

(5.12c) =0 on Se 

admits only the trivial solution -tjj = 0. 

The main difhculty in the proof of Lemma 15.51 is that the domain depends on 
the small amplitude parameter e > whereas the statement of Lemma 15.51 calls for 
an estimate of solutions of the system (I5.12p uniform for e > 0. In order to compare 
(I5.12p for different values of e, we employ the action-angle variables, which map the 
domain into a common domain independent of e. For any {x,y) € D^, let us 
denote by {{x',y') : ^'^ix' ,y') = ^'^i^TU) = p} the streamline containing {x,y), by 
a the arc-length variable on the streamline {V'e(a;',y') = p}, and by a{x,y) the 
value of a corresponding to the point (a;, y) along the streamline. Let us define the 
normalized action-angle variables as 

(5.13) 

1=—-^ dy'dx', and 9 = v,{I) —— da', 

J J {^,{x' ,v')<p} Jo V^Vt] {■>PAx'y')=p} 

where /iq and /if are the mean water depth at the parameter values and e, respec- 
tively, and 
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The action variable / represents the (normaUzed) area in the phase space under the 
streamhne {{x',y') : ipe{x',y') — '0e(a:,y) — p} and the angle variable 6 represents 
the position along the streamline of ipg{x,y). The assumption of no stagnation, 
i.e. ipeyix^y) < throughout V^, implies that all stream lines are non-closed. For, 
otherwise, the horizontal velocity ipey must change signs on a closed streamline. 
Moreover, for e > sufficiently small, all streamlines are close to those of the trivial 
flow, that is, they almost horizontal. Therefore, the action-angle variable {6, 1) is 
defined globally in V^. The mean-zero property (|2.8p of the wave profile r]^{x) 
implies that the area of the (steady) fluid region is {2TT/a)h^. Accordingly, by 
the definition in (|5.13[) it follows that < I < Hq and < 9 < 2TT/a, independently 
of e. 

Let us define the mapping by the action-angle variables by Ae{x,y) = (0,1). 
From the above discussions follows that Ae is bijective and maps to 

D ^ {{9,1) ■.0<9 < 27r/a, < / < ho}. 

At e = 0, the action-angle mapping reduces to the identity mapping on Vq ~ D. 
For e > 0, the mapping has a scaling effect. More precisely, 



j))d9dl = ^ II f{x,y)dydx 



ID J jj, 

for any function / defined in 

Another motivation to employ the action-angle variables comes from that they 
simplify the equation on the particle trajectorjQ. In the action-angle variables {9, /), 
the characteristic equation ()3.8|) becomes ([4, Section 50], [HI p. 94]) 

\/ = 0, 

where the dot above a variable denotes the differentiation in the (T-variable. This 
observation is very useful for future considerations. 

Proof of Lemma 15.51 Suppose on the contrary that there would exist sequences 
Efe 0+, Afc ^ Ao as fc ^ oo and V'fc G H^{T^€k) such that -0^ ^ is a solution of 
(|5.12p with e — ek- After normalization, ||V'fc||L2(r>,j,) = 1- We claim that 

(5.14) WMhhv^j^C, 

where C > is independent of k. Indeed, by Minkovski's inequality it follows that 



(5.15) 



= ll7'(^.JIlL^||Vfe||L^(p.j (i + ^1^) ^ II7'(V'.JIIl^ (i + M±i j . 

This uses the fact that the mapping (x, y) — > {X^^ (s), Y^^ (s)) is measure preserving 
and that |A — Ao| < ReAo/2. The standard elliptic regularity theory [2 for a 



■^In the irrotational setting, i.e., 7 = 0, a qualitative description of particle trajectories is 
obtained |16) by studying a specific nonlinear boundary value problem for harmonic functions. 
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Neumann problem adapted for (|5.12|) then proves the estimate (I5.14|) . This proves 
the claim. 

In order to study the convergence of {ijjk}, we perform the mapping by the 
action-angle variables (|5.13p and write Ae^{x,y) = {0,1). Note that the image of 
2?e^ under the mapping Ae^. is 

D = {{9J) -.0 < 6 < 2TT/a , < / < ho}, 

which is independent of k. Let us denote Aipk{0,I) = ipk{A'^^{d, I)). It is immedi- 
ate to see that A-ipk G H^{D). 

Since = ho + 0{e) it follows that 

WA^kWmo) - {ho/h,)Uk\\mv,j - 1 + 0(e). 
This, together with (I5.14p . implies that 

A^k "000 weakly in H^{D) as k ^ oo, 
A^k "000 strongly in L^{D) as fc — > oo 

for some 0oo- By continuity, ||'0oo|lL2(£)p') ~ 1. Our goal is to show that ipoo = 
and thus prove the assertion by contradiction. 

At the limit as — > 0, the limiting mapping is the identity mapping on 
D — Vq and the limit function ip^o satisfies 

A0oo+7'(0o)V'oo-7'(V'o) / Xoe^'>'iPooiXois),Yo{s))ds = in D; 

-/ — OO 

(5.16) dy'tlJoo{x,ho) = 0; 

000 (a;, 0) =0. 

Here, Aq = -iaca and {Xq{s),Yo{s)) = {x + U{y)s,y). 

Since the above equation and the boundary conditions are separable in the x 
and y variables, -000 can be written as 

OO 

V'oo(x,2/)=5]e*'"-0,(y). 

In case I — 0, the boundary value problem ()5.16|) reduces to 

0^,' = for ye (0, ho) 
(b'oiho) = 0, 0o(O) = 0, 

and thus, 0o = 0. 

Next, consider the solution 0i of (|5.16p when 1 = 1: 

0'/ - a20i - ^01 = for 2/e(0,/io) 
0i(/^o) = O, 0i(O) = O. 

This uses that j'itpo) = —U"/U. Recall that for the unstable wave number a and 
the unstable wave speed Cq, there exists an unstable solution 0a of the Rayleigh 
system (|5.4p . As is done in the proof of Lemma [5^ an integration by parts yields 
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that 



= / ( 01 ( 0a - «'0« - jt-TtJ-) - '^^ - "'-^i - [T^'^i ' ' '^^ 



0q(^o)0i(^o)- 



,(C/(M-c„)2 t/(/^o)-c„. 

Since 0q(/io) 0, we have 0i(/io) = 0; which together with 4i'i{ho) ~ 0, imphes 
that 01 = 0. 

For Z ^ 2, the solution 0; of (|5.16p ought to satisfy 

U[' -Pa^i- jj^^cf^i^O for ye(0,/io) 
\0;(/io) = O, 0i(O) = O. 

An integration by parts yields that 

ho , _ JJII 



and subsequently, for any q real it follows that (see the proof of Lemma [4?6l 

The same argument as in proving Lemma 14.61 applies to assert that 

;2„2 



Jo 

Recall that K{y) = -U"{y)/{U{y) - U,) > 0. 

Let —a^ be as in Lemma [5^ the lowest eigenvalue of — K{y) on y £ (0, ho) 
with the boundary conditions (f>'{ho) = and 0(0) = 0. By Lemma l5.41 amax > ctn- 
On the other hand, the variational characterization of — asserts that 

\\<P['-K{yMi\^)dy;^-al \<Pi\'dy. 
Jo Jo 

Accordingly, 

''''"{Pa'-alMil'dy^O 



must hold. Since a > Q!max/2 and Z ^ 2, it follows that Pa^ — an > 2a^-,^^ — > 0. 
Consequently, (f)i = 0. 

Therefore, V'oo = 0, which contradicts since H^ooHl^ = 1- This completes the 
proof. □ 

Lemma 5.6. Under the assumption of Lemma \5.5[ for any b E L'^{Se), there exists 
an unique solution ipi, to |5. ?| ). Moreover, the estimate 

(5-17) \\MHHv,)^C\\bU2^s^j 

holds, where C > is independent of e, b. 

Proof. The proof uses the theory of Fredholm alternative as adapted to usual elliptic 
problems [28l Section 6.2]. Let us introduce the Hilbert space 

(5.18) H{V,)^{ijeH\V,):iP{x,0)^0} 
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and a bilinear form Bz : H x H defined as 

Here, 

.0 

z e M and (•, •) denotes the (V^) inner product. By the estimate (|5.15p it follows 
that 

(5.19) s=; iC{j,S) + zm\H4^\\m 
and 

(5.20) ^ coll 
for z > 2C(7, 5), where 

and Co = min(l, C(7, (5)) > 0. Then, by the Lax-Milgram theorem there exists a 
bounded operator : H* H such that Bz[Lzf, (j>] = (/, <j>) for any / £ H* and 
(j> G H, where H* denotes the dual space of H and (•, •) is the duality pairing. For 
b G Lp'iSf), the trace theorem [351 Section 5.5] permits us to define b* e H* by 

<b*,il;>= <j) b{x)\p*{x,ri^{x))dx for -0 e iJ. 

Note that ipt, E H is a. weak solution of (|5.7p if and only if 

Bz{i!b,(l>) =< zipb + b* ,4)> for aU (/) e iJ. 
That is to say, ipb = Lz {zipb + b*), or equivalently 

(5.21) {id- zLz)i^b ^ Lzb*. 

The operator Lz : LF'{T>f) L'^{'Df) is compact. Indeed, 

\\LMm{v,) ^ ll-^zll//*^// II'/'IIh- < C'||</>|1^2(i,^) 

for any (/i G L?'{T>f). Moreover, the result of Lemma ISTSl states that ker(/ — zLz) = 
{0}. Thus, by the Fredholm alternative theory for compact operators, the equation 
(|5.2ip is uniquely solvable for any b e L'^{Se) and 

(5.22) i^b = {id- zLzr'^ Lzb*. 
Next is the proof of ([5T7)) . From it follows that 

Hbh^v,) ^ \\{-id- zLzy^\\^2^L2 \\Lz\\h'^h P*\\h' ^ C\\b\\L^(^s,), 
where C > is independent of b and e. Then, by (|5.2ip it follows that 

W'^bWH^V.) = \\Lz{ztpb) +Lzb*\\fj < C'dl^'bllL^CD,) + \\b\\L^(s,)) SS 

where C, C" > are independent of b and e. This completes the proof. □ 

Similar considerations to the above proves the unique solvability of the inhomo- 
geneous problem. 
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Corollary 5.7. Under the assumption of Lemma \5.5\. for any f G H*['Dt) there 
exists an unique weak solution ?/'/ G H^(T)c) to the elliptic problem 

Ai^ + 7'(V',)V-7'(V'.) / Xe^'^P{X,{s),Y,{s))ds^f mV,, 



ipn = on Se, 
tpix,0) = 0. 

The estimate 

(5.23) Uf\\m(v,)^C\\f\\H*(v,) 

holds, where C > is independent of f and e. Moreover, if f ^ L^iV^) then an 
improved estimate 

(5.24) Uf\\HHv._)^C\\f\\L^v,) 
holds, where C > is independent of f and e. 

Proof. As in the proof of Lemma 15. 6i the function ipf is a. solution to the above 
boundary value problem if and only if 

{id ~ zLz)^f = Lzf. 

The unique solvability and the i/^-estimate (|5.23p are identically the same as those 
in Lemma [5.61 When / e L^{'D^) the elliptic regularity theory 2 for Neumann 
boundary condition implies (|5.24p . □ 

For any b E Lp^^.{Se) let us define the operator 

(5.25) %b = ipbls, = ipbix,ri^{x)) , 

where ipb is the unique solution of (|5.7p in Lemma 15.61 with periodicity. Then, the 
elliptic estimate (|5.17p and the trace theorem PSI Section 5.5] implies that 

(5.26) \\m\m/2(s,) < CUbWmiv,) ^ C'Wbh^; 
Therefore, % : Lp(,,.(5e) Lp^j.{Se) is a compact operator. 

5.3. Proof of Theorem 15.11 This subsection is devoted to the proof of Theorem 
15.11 pertaining to the linear instability of small-amplitude periodic traveling waves 
over an unstable shear flow. 

For e ^ sufhciently small and |A — Ao| ^ (ReAo)/2, let us denote 

(5.27) T{X, e) = r (A, e) C (A, e) {P,yC (A, e) + n id) : Ll^^{Se) ^ i^er(5.), 

where C(A,e) = and T(A, e) = % are defined in (|5.5p and (|5.25p . respectively. 
In the light of the discussion in the previous subsection, it suffices to show that for 
each small parameter e > there exists A(e) with |A(e) — Ao| ^ (ReAo)/2 such that 
the operator id + JF(A(e), e) has a nontrivial null space. The result of Theorem 14.21 
states that there exists Aq = —iaCa with Imca > such that id + .F(Ao, 0) has a 
nontrivial null space. The proof for e > uses a perturbation argument, based on 
the following lemma due to Steinberg [STj . 

Lemma 5.8. Let F{X, e) be a family of compact operators on a Banach space, 
analytic in X in a region A in the complex plane and jointly continuous in (A, e) for 
each (A, e) e A x R. Suppose that id — F{Xo,e) is invertible for some Ao G A and 
all e e R. Then, R{X,e) = {id — F{X,e))^^ is meromorphic in A for each e G R 
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and jointly continuous at (zo,eo) Aq is not a pole of R{X,e); its poles depend 
continuously on e and can appear or disappear only at the boundary of A. 

In order to apply Lemma l5.8l to our situation, we need to transform the operator 
(|5.27p to one on a function space independent of the parameter e. This calls for the 
employment of the action-angle mapping At, as is done in the proof of Lemma l5.5l 

Ae'-V.^D and Aeix,y) ^ {9, 1), 

where the action-angle variables {9, 1) are defined in (|5.13|) and 

D = {{9, 1):0<9< 27r/a , < / < ho}. 

Note that Ae maps bijectively to {{9, ho) : < 9 < 27r/a}. The latter may be 
identified with (27r/a, ho). This naturally induces an homeomorphism 

B,:Ll^,{S,)^Ll^^{[0,2n/a]) 

by 

{BJ){9)^f{A-H0,ho)). 
Let us denote the following operators from L^^^. ([0, 27r/a]) to itself: 

f{X,e)^B,T{{X,e)){B,)~\ 
C{X,e)^B,C{X,e) {B,)-\ 
V, ^ B,P,y{B,)-\ 

and 

(5.28) T{X,e) ^B,T{X,e){B,)-^ =f {X,e)C{X,e){V,C {X,e) + VLid). 

Since T(A, e) is compact and T(A, e) and C(A, e) are analytic in A with |A — Ao| ^ 
(ReAo)/2, the operator !F{X^ e) is compact and analytic in A. Subsequently, !F{X^ e) 
is compact and analytic in A. Clearly, P^y{x), C(A, e) and B^ are continuous in e, 
and in turn, C(A, e) and are continuous in e. The key technical lemma is to show 
the continuity of T(A, e) in e and thus obtain the the continuity of !F{X, e) in e. 

Lemma 5.9. For e ^ sufficiently small and |A — Ao| ^ (RcAo)/2, the operator 
T(A, e) satisfies the estimate 

(5.29) ||T(A, ei) - T(A, e2)\\^^^xio.'2n/a])^^^j[o.27v/a]) ^ C\ei - 62!, 
where C > is independent of X and e. 

Proof For a given b £ ip„.([0, 27r/Q!]), let us denote b^ = B^^b e L^{Se). By 
definition 

f {X,e)b ^ BM = B.iiPbAsJ, 
where ipb.e & (^c) is the unique weak solution of 

(5.30a) Ai;b,e+l'{A)i'b,e-j'{iJ.) Xe^'MXe{s),Y,{s))ds = in P,; 

•/ —00 

(5.30b) {ipb,e)n = h on S^; 

(5.30c) Vf.,e(a;,0) = 0. 

Our goal is to estimate the i^-operator norm of T(A, ei) — T(A, 62) in terms 
of |ei — £2!. Since the domain of the boundary value problem (|5.30p depends on 
e, we use the action-angle mapping Ae^ {j = 1, 2) to transform functions and the 
Laplacian operator in (j = 1, 2) to those in the fixed domain D. To simplify 
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notations, we use Ae^ {j ~ 1, 2) to denote the induced transformations for functions 
and operators. Let ipj — Aej{ipb.ej), which are i/^(Z?)-functions, and let 

A,=A,(A), 7,(/)=A,(7'(^e,)), 

where j — 1,2. By definition, {b^ ) = B^ .{b^ .) — b. Note that the characteristic 
equation (|3.8p in the action-angle variables (9, 1) becomes 

\i = o 

for j = 1,2, where Vj{I) = v^.{I). That means, the trajectory {X^-{s)^Y^.{s)) in 
the phase space transforms under the mapping A^^ into {6 + Vj{I)s,I). Since 4'j 
(j = 1, 2) are 27r/a-periodic in 6, we have the Fourier expansions 

Under the action-angle mapping Ae- (j ~ 1, 2) the left side of (I5.30p becomes 



(5.31) =7,(/)(Ee^'V,./W" r Ae^^^ 



A + ilvj{I) 
and thus the system (|5.30p becomes 

dii>,{eM) = b{e), 

V',(0,O) = O. 

Accordingly, the difference -01 — V'2 is a weak solution of the partial differential 
equation 

(5.32) 

Ai(V'i - i^2) + (Ai - A2)^2 + 71 E T^^^iM^) - 

^ A + llVi 

+ 71 E (j^ - J^) Miy'' + (71 ~ 7.) E J^MDe^'' = 
with the boundary conditions 

(5.33a) 9/(01 - V2)(0, /^o) =0, 

(5.33b) (V^i-V2)(^,0) = 0. 

Let us write (|5.32p as 

(5.34) Ai(0i - 02) + 71 E A-^Ifci?(/) ^^^''^^^ ~ V'2,i(/))e'''' = /, 
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where / = /i + /2 + /s with 

/l =-(Ai-A2)V2, 

We estimate /i,/2,/3 separately. To simphfy notations, C > in the estimates 
below denotes a generic constant independent of e and A. 
First, we claim that /i e H*{D) with the estimates 

(5-35) ^ " e2|||fo||L2([0,2Va])> 

where H*{D) is the dual space of 

H{D) = {V'G : V(6l,/) = V(6' + 27r/a,/), V(6',0) =0}. 

Let us write 

Aj = a^jjdii + a'jgdie + a^gdee + b'idi + bld0 for j ^ 1,2, 
and the difference of the coefficients as 

- 12 - 12 - 12 

ail — 0/7 — Ojj, aie — a^g — a^g, age ~ agg — agg, 
bi^b]- bj, be^bl- bl 
Then formally, for any cj) £ H{D) D C^i^D) it follows that 

(5.36) ( M dide 

J Da 

-~4> {audu + ajedie + aggdee + bgdg + bjdi) ip2dld9 

Do 

[diijj2di (aiicj)) -I- diip2dg (aigcp) + dgip2dg {agg(l))] dIde 

Da 

<t> {bidii^2 + bgdg^P2) dldO ~ [ a//06 {0) dO, 

Da J{I=ho} 

This uses that and (j) are periodic in the ^-variable and that 

diM9,ho)=b{9), 0(0,O)=O. 

Note that the elliptic estimate (|5.17[) and the equivalence of norms under the trans- 
formation A^^^ assert that 

MhHD) ^^\\-^e2^M\m{V,,) = II V'b.e. 11^1(15, J 
IIl2(5,2) ^ ^ ll^llL2([o,27r/a]) ' 

Since 

\aii\ci + \aie\ci + l^eelci + \^e\ci + l^-f|ci ^ '^d^i ~ ^^l), 
by using the trace theorem it follow from (15.361) the estimate 



fi<j) dldO 

D 



< C |ei - 62! ^^2 II 7/1 + ll&llL2([o,2Va])) II'?^IIh1(D) 

< C |ei - £21 ll^llL2([o,2Va]) U\\m(D) ■ 
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This proves the estimate (|5.35p . 

We claim that if b is smooth then the formal manipulations in ()5.36|) are valid and 
■02 S C'^{D). Note that Theorem l2 . 2l ensures that the steady state [rj^^ (x), ij]^^ {x, y)) 
IS m C^+P class, where /3 G (0, 1). Since 6 is smooth it follows that b^^ = b 
is at least in (^ej). Then, the similar argument as in the regularity proof of 
Theorem O below asserts that Vb.ea & ^^^^^(^'es) C C^iV^^). Since the definition 
of the action- angle variables guarantees that the mapping is at least of C^, 
subsequently, '!/'2 = -^t^'^b.ei ^ C^{D). This proves the claim. If 6 e L^, an 
approximation of b by smooth functions establishes (|5.35p . 

Next, since 



1 



A -I- ilvj 



1 1 



(IReAP + IImA-^Wjf )i/2 ^ |ReA| " ReAo' 
by the estimate (|5.17p it follows that 

(5-37) ||/2||l2(d) 5^ C\ei - e2|||V'2|lL2(D) < C\ei - e2|||6||L2([o,2,r/Q])- 

Similarly, 

(5-38) Wfah^Do) ^ '^l^i - <^2\\\b\\L^[a,2^/a])- 

Combining the estimates (I5.35p . (|5.37p and (|5.38p asserts that / £ H*{D) and 

\\f\\H*(D) ^ C'kl - e2|||fe||L2([0,277/Q])- 

Let ^ = "01 — "02 G H^{D) and = A'^^^tp G (T^ei)- It remains to transform 
back to the physical space of the boundary value problem for and to compute the 
operator norm of T (A, ei) — T (A, £2). Under the transformation Aj^ , the equations 
(|533a|) . (|5.33bp become 



J —00 

(0^On=O on 5,,; 

Then, A-^f e H*{V,^) and 

l|-^iriV||jj.(27^j ^ C \\f\\H'{D) ^ C:\ei - e2|||6||L2([o,2,r/Q])- 

Corollary [577] thus applies to assert that 

MmiD) < M\m{v,,) < C II A"; V II (P. J 

< C|ei - e2|||&||L2([o,27r/a])- 

Finally, by the trace theorem it follows 
||T(A,ei)6-r (A,e2)6|U2^^([o,2V"]) = IKV'i " ^2) ^o)lli2^^([o,2Va]) 

< \\tp\\Hi{D) ^ C'kl ^ e2|||^'||L2^„,([0,27r/Q])- 

This completes the proof. □ 
We are now in a position to prove our main theorem. 
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Proof of Theorem \5.1\ For |A — Ao| ^ (ReAo)/2, where Aq ~ —iaca, and e ^ 
small, consider the family of operators ^{X, e) on Lp^^{[0, 2n/a]), defined by (|5.28p . 
The discussions following (|5.28p and Lemma 15.91 assert that !F{X, e) is compact, 
analytic in A and continuous in e. 

By our assumption, Imca > and Ca is an unstable eigenvalue of the Rayleigh 
system (|4.1l) - (|4.2p which corresponds to e = 0. In other words, Aq is a pole of 
{id + T{X,0))~^ . Subsequently, it is a pole of {id + !F{X,0))~^ . Since Aq is an 
isolated pole, we may choose S > small enough so that the operator id + J-{X, 0) 
is invertible on |A — Ao| = (5. By the continuity of J-{X, e) in e, the following estimate 

\\T{X,e) - ^(A,0)||i2^_^([o^27r/Q])-*L|„([0,27r/a]) ^ Ce 

holds. Hence, id + J-{X, e) is invertible on |A — Ao| = 5 for e ^ sufficiently small. 
Then by Lemma l5.8[ the poles of {id-\- T{X, e))^^ are continuous in e and can only 
appear or disappear in the boundary of {A : | A— Aq < (5}. Therefore, for each e ^ 0, 

there exists a pole A(e) of {^d + T{X,e)^ in |A(e) — Ao| < (5. Thus, ReA(e) > 

and there exists a nonzero function / e Lpj,j.([0, 27r/a]) such that 

(5.39) (id + J^(A,e))/ = 0. 

Below we construct an exponentially growing solution to the linearized system 
(f3:6a|l - ([3Je| . Define 

f^{B,r'feL\S,), 

then 

(5.40) (/ + .F(A,e))/-0 

by ()5.39p and the definition of T. Let -0 {x, y) E (I?e) to be the unique solution 
of dnZl) with A = A(e) and 

(5.41) i^n{x) = b = -C\P,y{x)C^ + ni)f{x). 
By ([00]) . ((5^ and the definition of we have 

f ^%b^ ilj{x,f]^{x)) 

and thus 

0„(X) = -C\P,y{x)C^+ni)^^{x,7],{x)). 

Define 

(5.42) r,{x) = mx, r,,{x))] G L^{S,), 
and 

P{x,%{x)) = ~Pey{x)r]{x), 

then (|5.4ip becomes 

(5.43) = -C^ (Pi^, Vc{x)) + niP{x, rj,{x))) . 

Now we show that [e'^^'^'^* ijj {x, y) , e^^'^^*^r]{x)\ satisfies the linearized system p.6ap - 
p.6el) . The bottom boundary condition (|3.6ep is satisfied since ip{x,0) = 0. The 
equation (|3.6cp is automatic. By Lemma [531 and equations (|5.42l) . (I5.43p . t] {x) and 
ipn{x) satisfy the equations (|5.1bp and (j5.1dp weakly. Equivalently, the equations 
p.6bp and (|3.6dp are satisfied weakly. Since 'tp{x,y) satisfies the equation (|5.7b ). 
we have 
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(5.44) w = -AV' = 7'(V'.)V' - 7'(V'.) / Ae^XX,(s),y,(s))ds. 

As shown in |42| , above equation implies that the vorticity lu satisfies the equation 
(|3.7p weakly. Equivalently, the equation (j3.6ap is satisfied weakly. In summary, 
|^gA(e)t^ (x, y) , e'*'^'^*77(a;)] is a weak solution of the linearized system (|3.6a[) - (|3.6ep . 

Our last step of the proof is to get the regularity of the growing-mode {e^'-'^'>*'r]{x), 
e^^'^^*^{x, y)) and thus show that it is a classical solution of (|3.6a[) - (|3.6ep . By (|5.4ip 
and Lemma Em it follows that tp e H'^ (Pg). We claim that ^ e {V^). Indeed, 
by the trace theorem, ip e (V^) implies that ip{x,ri^{x)) £ iJ' (iS^). Since the 
operator is regularity preserving, by (|5.43p tpnix) £ [Sg). This, together 
with the facts that lo ~ — AV' £ {T^e) and that the steady state 

(77,(0:), V. (a:,2/))eC3+",ae (0,1) 

(see [IH] or Theorem [^^ . implies that ip £ (V^) by the regularity theory ([5]) 
of elliptic boundary problems. Then, by using the trace theorem and (|5.43p again, 
we get ip{x,ri^{x)) and ?A„(a;) £ [S^). 

In order to obtain the higher regularity for ■(/;, we need to show that lo £ H^{T>^). 
The argument presented below is a simpler version of that in [40) . Taking the 
gradient of (|5.44p yields that 

Vc^ =V(7'(Ve))V' + 7'(V'.)VV' - V(7'(V'e)) / Ae^XX,(s), y,(s))ds 

(5.45) „ -'-^ 

- 7'(^.) / Ae^^V^(X,(.), y,(.))^i^^M4Md,. 

Note that the particle trajectory is written in the action-angle variables (0, /) = 
Ae{x,y) as 

(X,(s; X, y), y,(s; x, y)) ^ A-;\(6 + v,{I)s, I)). 

This relies on that the action-angle mapping Ae is globally defined, as a consequence 
of the fact that the steady flow has no stagnation. With the use of the above 
description of the trajectory the estimate of the Jacobi matrix 



(5.46) 



d{X,{s;x,y),Y^{s;x,y)) 



< Ci \s\ + C2 



d{x,y) 

follows, where Ci, C2 > are independent of s. 

It is straightforward to see that by calculations as in proving (|5.15p . the L^-norm 
of the first three terms of (|5.45p is bounded by the norm of t/j. The last term 
in (|5.45p is treated as 



7'(^.)/° Ae^^V^(X,(.),y,(.))&^)4M)d, 



L2 



d{x,y) 

^ ll7'(V'e)||L^ / \X\e^^^'{C^\s\+C2m^/J{X,{s),Y,{s))\\L2^^^)ds 



UNSTABLE SURFACE WAVES IN RUNNING WATER 



53 



This uses (|5.46p . ReA ^ S > and the fact that the mapping (x, y) ^ (Xe(s), i^e(s)) 
is measure-preserving. Therefore, 

In turn, uj e H^{'D^). Since ipnix) <E H^^^{Se)^ by the elhptic regularity theorem 
[2j it follows that that e {T^e)- In view of the trace theorem this implies 

We repeat the process again. Taking the gradient of ()5.45|) and using the linear 
stretching property (??) of the trajectory, it follows that uj S [T^t)- The elliptic 
regularity applies to assert that ip G H^{V^) C C'^+^{'D^), where (3 € (0, 1). By the 
trace theorem then it follows that 'ip{x, r]^{x)) € H''/'^{S^). On account of (|5.42|) this 
implies that 77 e H^^ (S,) C C^+>^{[0,2Tr/a]). Therefore, (e^(')*77(x), e^(^)*i/'(a;, y)) 
is a classical solution of (|3.6[) . This completes the proof. □ 

6. Instability of general shear flows 

Linear instability of free-surface shear flows is of independent interests. This 
section extends our instability result in Theorem 14.21 to a more general class of 
shear flows. The following class of flows was introduced in [39] and [42] in the 
rigid-wall setting. 

Definition 6.1. A function U E C'^{[0,h]) is said to be in the class T if U" takes 
the same sign at all points such that U{y) — c, where c is in the range of U but not 
an inflection value ofU. 

Examples of the class- flows include all monotone flows and symmetric flows 
with a monotone half. Moreover, if U"{y) — f{U{y))k{y) for / continuous and 
k{y) > 0, then U is in class All flows in class IC~^ are in class J-. 

The lemma below shows that for a flow in class J- a neutral limiting wave speed 
must be an inflection value. The main difference of the proof from that in the class 
JC~^ case (Proposition l4.4|) is the lack of an uniform i/^-bound for the unstable mode 
sequence. 

Lemma 6.2. For U E J- , let {((/'fc, afe, Cfe)}^]^ with Imcfc > Q be a sequence of 
unstable solutions satisfying j^. If{ak,Ck) converges to {as,Cs) as k ^ 00 

with as > and Cg is in the range of U , then Cg must be an inflection value of U . 

Proof. Suppose on the contrary that Cg is not an inflection value. Let j/i, j/2, ■ • ■ , J/m 
be in the pre- image of Cg so that U{yj) = Cg, and let Sq be the complement of 
the set of points {yi, y2, ■ ■ ■ , ym} in the interval [0, h]. Since Cg is not an inflection 
value. Definition 16.11 asserts that U"{yj) takes the same sign for j — 1,2, ... ,m, 
say positive. As in the proof of Proposition 14.41 let Eg = {y E [0, h] : \y — yj\ < 
5 for some j, where j = 1, 2, • • • , m}. It is readily seen that Eg C 5*0. Note that 
U"{y) > ior y E Eg a S > small enough. We normalize the sequence by setting 
||'^fc||L2 = 1- The result of Lemma 3.6 in ,39 J implies that 0^ converges uniformly 
to (f>g on any compact subset of Sq. Moreover, 0" exists on So and (pg satisfies 

~ a^gPg - j^^cPg = for ye{0,h). 
U - Cg 

Our first task is to show that (pg is not identically zero. Suppose otherwise. The 
proof is again divided into two cases. 
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Case 1: U{h) ^ Cg. In this case, [h — 5i, h] C 5o for some 5i > 0. As is done in 
the proof of Proposition 14.41 for any q real, it follows that 

l-h 



h 



Jo J El Cfcl J^^ \U — Ck\ 

^ rKrdy + a^-sup '^7^~f / m^dy. 

Jo El \U~Ck\'' J El 

We choose q = J/min - 1, then by (|4.10p 
(6.1) £ (l^'.l^ + a?-|0.P + 

Imcfe / 



If e is chosen to be small then the above two inequalities lead to 
(6.2) 0^«g-sup '^"^^~^-'" + ^^' / \cb,\'dy-cj' \^,\'dy. 



\U-Ck\' 



Since (j)k converges to = uniformly on Eg and [h — di,h], this implies ^ as/2 
when k is large enough. A contradiction proves that (f>s is not identically zero. 
Case 2: U{h) = Cg. From ()4.2ip . we have 

" + 2al + at dy = 2alRegM) ih)\' + f i^^J^llM^dy. 

^ ^ Jo \U - Ck\ 

The imaginary part of (|4.16p yields 

(6.3) Imcfc ( y ^^^\2 dy = -Img,j(cfe) {h)\^ . 

Jo \U~Ck\ 

Denote C/^'ax = niax[Q /jj U"{y). Combining the above two identities, we have 



(6.4) im'^2alm' + atW?)dy^ / ^ ^^-,7 ^ dy 

Jo ^ 'Jo \U - Ck\ 

= {2alReg.{cu) - [U'^^. + 1)) (/.)|^ 

^ C'd (cfe, C/(/i)) (/^)|' Cd (cfe, C/(/i)) , 
where we use (|4.20p . Since 

dick,U{h)) ^ iRecfe - V{h)\ + (Imcfc)2 0, 
so for fc large enough we have 



2 i^l 
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which is a contradiction. This proves that (ps is not identically zero. Subsequently, 
Lemma 14.71 asserts that 0s [Vj ) 7^ for some j/j . 

Below, we get a contradiction from the assumption that Cg is not an inflection 
value. In Case 1 when U{h) ^ Cs, it is straightforward to see that 

'-^ ' dy^ / — -\(l3s\ rfy = oo, 



since <f>s{yj) 7^ and by our assumption U" > on {\y — yj\ < S}. Fatou's lemma 
then states that 

hmmf / — \(f)k\ dy = 00. 

Then similar to the estimate (16.11) above, we have 

n f'^ , 2 1. 12 , U"{U-Unun + l) ,^ |2^ , 

Jq V +"'=''^'=1 + \u ~-c\^ ) ^ 

Regrick) + (Rec, - C/,„in + |0,(/,)|2 

Imcfc / 

^ 177 1^ M'dy-sup — Ce>0, 

for k large. A contradiction asserts that Cs is an inflection value. For Case 2 when 
U{h) = Cs, similarly we have 

hm mf / — ^2 |0fc| dy^+oo, 

Jes \U-Ck\ 

and from (16.41) 



0^/ ^ dy-sup ^ >0, 

Jes \U-Ck\ El \U-Ck\ 

when k is large. Another contradiction completes the proof. □ 

The proof of above lemma indicates that a flow in class J- is linearly stable when 
the wave number is large. 

Lemma 6.3. Assume g =/= 1. Then for any flow U (y) in class there exists 
Q!max > such that when a ^ amax there is no unstable solutions to Il4.1\l~(4-Ml- 



Proof. Suppose otherwise. Then, there would exist a sequence of unstable solutions 
{{4>k, Oik,Ck)}^^i of (|4.ip - (|4.2p such that a/c ^ cx) as fc — > oo. After normaHzation, 
let ||0fc||L2 = 1- First we show that limfe^oo Imc/t = 0. If Imcfc ^ S > for some S, 
then I / \ U — Ck\ and |gr(cfc)| are uniformly bounded. Accordingly, with q = Reck 
in (|iT^ it follows that 

= (10^1' + al + I'^fcH dy - Rcgrick) \^k (Ml' 

Jo \ \U~Ck\ J 

^a^-sup ^^_^^^2 + l^l'kl dy-C[el \(j)k\ dy + - j |0fc| dy 

2 \U"(U~Reck)\ C „ 

\U~Ck\ £ 
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when k is big enough. This contradiction shows that Ck Cs (z [Umin, CAnax] when 
k oo. The remainder of the proof is nearly identical to that of Lemma 16.21 and 



The following theorem gives a necessary condition for the free surface instability 
that the flow profile should have an inflection point, which generalize the classical 
result of Lord Rayleigh 52J in the rigid wall case. 

Theorem 6.4. A shear flow U (y) without an inflection point is linearly stable in 
the free surface setting. 

Proof. Suppose otherwise; Then, there would exist an unstable solution a, c) to 
(|4.ip - (|4.2[) with a > and Imc > 0. Lemma l4!8l allows us to continue this unstable 
mode for wave numbers to the right of a until the growth rate becomes zero. Note 
that a flow without an inflection point is trivially in class F. So by Lemma 16.31 
this continuation must end at a finite wave number Ofmax and a neutral limiting 
mode therein. On the other hand. Lemma 16.21 asserts that the neutral limiting 
wave speed Cs corresponding to this neutral limiting mode must be an inflection 
value. A contradiction proves the assertion. □ 

Remark 6.5. Our proof of the above no-inflection stability theorem is very different 
from the rigid wall case. In the rigid-wall setting, where (j>{h) — 0, the identity 
(|4.15p reduces to 



which immediately shows that if U is unstable (c^ > 0) then U"{y) = at some 
point y € (0, /i). The same argument was adapted in [62i, Section 5] for the free- 
surface setting, however, it does not give linear stability for general flows with no 
inflection points. More specifically, in the free-surface setting, (I4.15P becomes 



which only implies linear stability ( |62[ Section 5]) for special flows satisfying 
U" (y) < 0, U' (y) ^ or U" (y) > 0, U' (y) ^ 0. In the proof of Theorem EM 
we use the characterization of neutral limiting modes and remove above additional 
assumptions. 

Let us now consider a shear flow [/ G with multiple inflection values Ui,U2, . ■ . ,Uj 
Lemma 16.21 states that a neutral limiting wave speed Cs must be one of the inflec- 
tion values Ui,U2, ■ . ■ ,Un, say Cg — Uj. By localizing the estimates in the proof 
of Lemma 14.61 around inflection points with the inflection values C/^, we can get an 
uniform bound for the unstable mode sequence. We skip the details, which are 
similar to the case of rigid walls treated in [42]. Thus, neutral limiting modes for 
flows in class J- are also characterized by inflection values. 

Proposition 6.6. IfU^T has inflection values Ui, U2, . ■ . , Un, then for a neutral 
limiting mode {(j)s, cts, Cs) with > the neutral limiting wave speed must be one 
of the inflection values, that is, Cg = Uj for some j . Moreover, (ps must solve 



hence is omitted. 



□ 





ct>'^ - alc^s + K,{y)c^s = for y ^ (0, h), 
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with boundary conditions 

where K,{y) - -U"{y)/{U{y) ~U,). 

One may exploit the instability analysis of Theorem 14.91 for a flow in class T 
with possibly multiple inflection values. The main difference of the analysis in class 
J- from that in class /C"*" is that unstable wave numbers in class T may bifurcate to 
the left and to the right of a neutral limiting wave number, whereas unstable wave 
numbers in class JC^ bifurcate only to the left of a neutral limiting wave number. 
In the rigid- wall setting, with an extension of the proof of Theorem 1.1], Lin 
[42l Theorem 2.7] analyzed this more complicated structure of the set of unstable 
wave numbers. The remainder of this section establishes an analogous result in the 
free-surface setting. 

In order to study the structure of unstable wave numbers in class T with possibly 
multiple inflection values, we need several notations to describe. A flow U ^ T \s 
said to be in class if each Kj{y) = —U"{y)/{U{y) — Uj) is nonzero, where Uj 
for j — 1, • • • , n are inflection values oi U. It is readily seen that for such a flow 
Kj takes the same sign at all inflection points of Uj. A neutral limiting mode 
{4>j, aj, Uj) is said to be positive if the sign of Kj is positive at inflection points of 
Uj, and negative if the sign of Kj is negative. Proposition lG. 61 asserts that —aj. is a 
negative eigenvalue of —j^ — Kj{y) on y e (0, h) with boundary conditions (|6.5p . 
We employ the argument in the proof of Theorem 14.91 to conclude that an unstable 
solution exists near a positive (negative) neutral limiting mode if and only if the 
perturbed wave number is slightly to the left (right) of the neutral limiting wave 
number. Thus, the structure of the set of unstable wave numbers with multiple 
inflection values is more intricate. We remark that a class-/C+ flow has a unique 
positive neutral limiting mode and hence unstable solutions bifurcate to the left of 
a neutral limiting wave number. 

Let us list all neutral limiting wave numbers in the increasing order. If the 
sequence contains more than one successive negative neutral limiting wave numbers, 
then we pick the smallest (and discard others) . If the sequence contains more than 
one successive positive neutral limiting wave numbers, then we pick the largest 
(and discard others). If the smallest member in this sequence is a positive neutral 
limiting wave number, then we add zero into the sequence. Thus, we obtain a new 
sequence of neutral limiting wave numbers. Let us denote the resulting sequence by 
<«([<•••< < a^, where, (might be 0),. . . , are negative neutral 
limiting wave numbers and aj, . . . , are positive neutral limiting wave numbers. 
The largest member of the sequence must be a positive neutral wave number since 
no unstable modes exist to its right. 

Theorem 6.7. For U G with inflection values Ui,U2, ■ ■ ■ ,Un, let < Qq < 
■ ■ ■ < < be defined as above. For each a G uj^Q(Q!j, a^), there exists 
an unstable solution of J^. JJ )~ |T^ Moreover, the flow is linear stable if either 
a ^ or all operators — Kj{y) (j = 1,2, . . . ,n) on y £ (0, h) with I16.5\ are 

nonnegative. 
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Theorem 16.71 indicates that there might exist a gap in (0, a^) of stable wave 
numbers. Indeed, in the rigid-waU setting, a numerical computation [7] demon- 
strates that for a certain shear-flow profile the onset of the unstable wave numbers 
is away from zero, that is > 0. 

Appendix A. Proofs of (fOSl) . (fOBl) . I[4AT\i and I^AA2h 
Our first task is to show that the limit (|4.35p holds as e 0— uniformly in 

In the proof of Theorem l4.9( we have already established that both (y; e, c) and 
(f>o {y; e, c) uniformly converge to 0^ in C^, as (e, c) (0, 0) in E(^ji bi.b2)- Moreover, 
(j)2{0;e,c) — > —-^qjyj uniformly as (e, c) (0,0) in -E(_R.fci.b2). So the function 

G{y, 0; e, c)(/)o(y; e, c) = (i^i {y; e, c) 02 (0; e, c) - (/)2 (y; e, c) c^i (0; e, c)) 0o(?/; e, c) 

convergesuniformly to— 03 (y) /^^ (0) inC^[0, d],as (e,c) ^ (0,0) in i?(fl_fcj^b2).Then 
the uniform convergence of (|4.35|) follows from (|4.33p . 
The proof of (j4.36p uses the following lemma. 

Lemma A.l ( 39 , Lemma 7.3). Assume that a sequence of differentiable functions 
{rfcjfe^i converges to Too 'in and that {ck}^i converges to zero, where Imck > 
and |Recfc| ^ Rlmck for some R > 0. Then, 
(A.l) 

U" K{y) ^ ^ Kjaj) ^ , 



- Jo [U-Us-Ckf " 7o tt|C^'(aj)l 

provided that U'{y) ^ at each Oj. Here ai, . . . ,am,, are roots of U — Us- 

We now prove (|4.36p . That is. We shall show that (I4.36P holds uniformly in 
^iRMM)- Suppose for some 5 > and a sequence {(£fc, Cfc)}^-^ in £'(fl f,i,()2)^it^^ 
max(6j,6|) tending to zero 

\ — {ek,Ck)-{C + iD)\>5 

holds, where C and D are defined in (|4.37p . Let us write 
(9$ /■'' U" d 

-^{£k,Ck) ^ - I 777 — 77 T^rfedy -t- — 5r(t^s + c)</'2(0;efc,Cfe) = /-H //, 

oc Jo [U -Us - CkY do 

where 

^k{y) = G{y,0;ek,ck)(l)o{y;ek,ck) --n-F7^4>l in ■ 

VsW 

By Lemma [A. 11 it follows that 

A straightforward calculation yields that 

hia II ^ -—gr[Us ' - ' ' - 



k^oo dc""' " 'cl)'M \{U{h)-Us)^ ({/(/i) - C/,)V 0^(0) 0',(O)' 

where A is given in (|4.3ip . Therefore, 

lim -^{ek,Ck) =C + iD. 

fc— >oo oc 
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A contradiction then proves the uniform convergence. 
The proofs of ()4.4ip and (I4.42|) use the foUowing lemma. 

Lemma A.2 ([39J, Lemma 7.1). Assume that {i/jk}'^i converges toipoo mC^([0,/i] 
and that {c^-}^]^ with Imcfe > converges to zero. Let us denote Wk{y) — 
U{y) — Us — Recfe. Then, the limits 



/loZd provided that U'{y) ^ at each aj. Here, ai, . . . , satisfy U{aj) — Us- 

We now prove gH]). Since K{y)(j>s(j>k -> K{y)(j)l in C^{[0,h]), by and 
(|X3|) it foUows that 



{U~Ck){U-Us) 



Ok 



4y 



(^•4) TT72^^^K{y)ct>sciokdy + i I -J^^^^K{y)<j,sAdy 

as fc ^ oo. Since Ck ^ Us as k oo in the proof of Theorem 14.21 it follows that 
(A.5) linr = .g;(C/.) = ^5 , C/'(/.) 



Addition of (|A.4p and (|A.5|1 proves (|4.4ip . In case U{h) — Us the same computa- 
tions as above prove (|4.42|) . This uses that 

lim g^(^fc) = _ ^(fe) - Cfc ^ p 

k^oo Ck — Us k^oo g + U'ih){U{h) - Ck) 
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